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Abstract 



In this article and the companion paper we address the question of how one might 
' obtain the semiclassical hmit of ordinary matter quantum fields (QFT) propagating on curved 

spacetimes (CST) from full fledged Quantum General Relativity (QGR), starting from first 
, principles. We stress that we do not claim to have a satisfactory answer to this question, rather 

our intention is to ignite a discussion by displaying the problems that have to be solved when 
carrying out such a program. 
^3 ■ III ths first paper of this series of two we propose a general scheme of logical steps that 

Q I one has to take in order to arrive at such a limit. We discuss the technical and conceptual 

O^' problems that arise in doing so and how they can be solved in principle. As to be expected, 

. completely new issues arise due to the fact that QGR is a background independent theory. 

For instance, fundamentally the notion of a photon involves not only the Maxwell quantum 
field but also the metric operator - in a sense, there is no photon vacuum state but a "photon 
^ ■ vacuum operator" ! Such problems have, to the best of our knowledge, not been discussed in the 

literature before, we are facing squarely one aspect of the deep conceptual difference between a 
background dependent and a background free theory. 

While in this first paper we focus on conceptual and abstract aspects, for instance the defini- 
tion of (fundamental) n— particle states (e.g. photons), in the second paper we perform detailed 
calculations including, among other things, coherent state expectation values and propagation 
on random lattices. These calculations serve as an illustration of how far one can get with 
present mathematical techniques. Although they result in detailed predictions for the size of 
first quantum corrections such as the 7-ray burst effect, these predictions should not be taken 
too seriously because a) the calculations are carried out at the kinematical level only and b) 
while we can classify the amount of freedom in our constructions, the analysis of the physical 
significance of possible choices has just begun. 
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1 Introduction 



Canonical, non-perturbative Quantum General Relativity (QGR) has by now reached the status of a 
serious candidate for a quantum theory of the gravitational field: First of all, the formulation of the 
theory is mathematically rigorous. Although there are no further inputs other than the fundamental 
principles of four-dimensional, Lorentzian General Relativity and quantum theory, the theory pre- 
dicts that there is a built in fundamental discreteness at Planck scale distances and therefore an UV 
cut-off precisely due to its diffeomorphism invariance (background independence). Next, while most 
of the results have so far been obtained using the canonical operator language, also a path integral 
formulation ("spin foams") is currently constructed. Furthermore, as a first physical application, 
a rigorous, microscopical derivation of the Bekenstein-Hawking entropy - area law has been estab- 
lished. The reader interested in all the technical details of QGR and its present status is referred to 
the exhaustive review article and references therein, and to for a less technical overview. For 
a comparison with other approaches to quantum gravity see |^, |^ . 

A topic that has recently attracted much attention is to explore the regime of QGR where the 
quantized gravitational field behaves "almost classical", i.e. approximately like a given classical 
solution to the field equations. Only if such a regime exists, one can really claim that QGR is 
a viable candidate theory for quantum gravity. Consequently, efforts have been made to identify 
so called semiclassical states in the Hilbert space of QGR, states that reproduce a given classical 
geometry in terms of their expectation values and in which the quantum mechanical fluctuations 
are small |0, |10|, |ll|. Also, it has been investigated how gravitons emerge as carriers of the 



gravitational interaction in the semiclassical regime of the theory |jT2|, |T3|, The recent investigation 



of Madhavan and others |T^, |T6|, |T^ on the relation between the Fock representations used in 
conventional quantum field theories and the one in QGR further illuminate the relation between 
QGR and a perturbative treatment based on gravitons. 

In this and the companion paper [Q] we would also like to contribute to the understanding of the 
semiclassical limit of QGR: We will investigate how the theory of quantum matter fields propagating 
in a fixed classical background geometry (QFT on CST) arises as an approximation to the full theory 
of QGR coupled to (quantum) matter fields. We will show in section ^ of the present work, how, 
upon choosing a semiclassical state, an effective QFT for the matter fields can be obtained from 
a more fundamental theory of QGR coupled to matter. This effective theory turns out to be very 
similar to standard QFT on CST, but still carries an imprint of the discreteness of the geometry in 
QGR as well as of the quantum fluctuations in the gravitational field. 

Validating the semiclassical limit of matter coupled to QGR is not the only motivation for the 
present work. Since QGR is a background independent theory, the consistent coupling of matter 
fields requires a quantum field theoretical description of these fields that differs considerably from 
that used in ordinary QFT. Therefore, another aim of the present work is to gain some insights 
into what these differences are and how matter QFT can be formulated in a setting where also the 
gravitational field is quantized. As a main result of the present paper we show how a theory of matter 
coupled to quantum geometry can be formulated within the framework of QGR. Within this theory 
we identify states that can roughly be compared to the n-particle states occurring in ordinary QFT. 
Their structure is however fundamentally different as compared to that of the ordinary Fock states: 
Their definition also involves operators of the gravitational sector of the theory! 

Similar considerations may be applied to understand the emergence of gravitons in the semi- 



2 



classical limit of QGR. The situation there is however a bit more complicated since it requires the 
separation of the gravitational field in a background- and a graviton-part. We refer the reader to 



19| where a detailed consideration will be given. 

Finally, due to a better understanding of the phenomenology of quantum gravity and the ex- 
periments that could lead to its detection (see pO[ for a recent review) it is an intriguing question 
whether it might already be possible to make predictions for observable quantum gravity effects 
based on QGR. In order to do so, one has to consider a coupling of the gravitational field to matter 
- one can not measure the gravitational field directly but only through its action on other fields. 
Indeed, 



2|. In 



ground-breaking work on the phenomenology of QGR has been done [ 
these works, corrections to the standard dispersion relations for matter fields due to QGR have been 
obtained. Since we are dealing with a theory for matter coupled to QGR in the present work, it is 



an important question whether the results of |21 



|, 3 can be confirmed in the present setting. 
We will discuss the general aspects of this question in section |^. In the companion paper [Q] we will 
carry out a more detailed calculation, based on the results of the present work and the semiclassical 
states constructed in IP, O, IIT 



The main difficulty in carrying out the program outlined up to now lies in the fact that the full 
dynamics of quantum gravity coupled to quantum matter is highly complicated. This would already 
be the case for ordinary interacting fields but is amplified in the present case due to the complicated 
interaction terms (the gravitational field enters in a non-polynomial way) and the difficulties in the 
interpretation of the resulting solutions. In the setting of QGR, the dynamics is implemented in 
the spirit of Dirac, by turning the Hamilton constraint of the classical theory into an operator and 
restricting attention to (generalized) states in its kernel. A mathematically well-defined candidate 
Hamiltonian constraint operator has been proposed in ^ |2^, ^ (see also ^ for 
another proposal based on Vasiliev invariants). This operator turns out to be very complicated and 
a systematic analysis of its kernel seems presently out of reach. Therefore, in our considerations, 
we can not start from a fully quantized dynamical theory of gravity coupled to matter. Instead, we 
have to treat the dynamics in some rather crude approximation and therefore our considerations 
will be kinematical to a large extent. To be more precise, we will not treat the the matter 
parts in the Hamiltonian as constraints, but as Hamiltonians generating the dynamics of the matter 
fields in the ordinary QFT sense. With the part in the Hamiltonian describing the self interaction of 
the gravitational field we will deal by using semiclassical states, which, as we will explain, annihilate 
this part of the Hamiltonian constraint at least approximately. Proceeding in this way certainly only 
amounts to establishing an approximation to the full theory: The self interaction of gravity and the 
back-reaction of the matter fields on the geometry are only partly reflected by using semiclassical 
states that approximate a classical solution to the field equations of the gravity-matter system. 

What we gain is a relatively easy to interprete, fully quantized theory of gravity and matter fields. 
This way we have "a foot in the door" to the fascinating topic of interaction between quantum matter 
and quantum gravity and can start to discuss the conceptual issues arising, as well as take some steps 
towards the prediction of observable effects resulting from this interplay. 

Let us finish this introduction with a brief description of the structure of the rest of the article: 
In the next section, we will discuss the main steps taken in this work in more detail before we turn 
to their technical implementation in the subsequent sections. 

In section ^ we give a very brief introduction to the formalism of QGR, mainly to fix our notation. 
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In more detail we display the matter Hamiltonian operators of electromagnetic, scalar and Dirac 
matter when coupled to general relativity. 

Section ^ contains the main results of this paper, namely a proposal for how to arrive at the 
notion of Fock states or n— particle states on fluctuating quantum spacetimes, if one is to start from 
a fundamental quantum theory of gravity of matter. 

In section ^ we discuss various methods to obtain dispersion relations for the matter fields from 
the full theory described in |^. 

We conclude this work with a discussion of its results and possible directions for future research 
in section |^. 

In an appendix, we treat the toy model of two coupled Harmonic oscillators to give an example 
of how the results are affected when one uses kinematical coherent states instead of coherent states 
in the dynamical Hilbert space of the theory. 

As already said, in the present paper we focus on describing the general scheme, detailed calcu- 
lations will appear in the companion paper 



2 A General Scheme 

In this section we want to discuss the issues related to the QFT on CST Limit of QGR and describe 
the steps taken in the present work in more detail. 



The first step that we will take is the kinematical quantization of the matter and the gravitational 
field on a Hilbert space T-^"^™. We will be guided by the fundamental principles of QGR which have 
to be obeyed: The quantum theory should be formulated in a background free and diffeomorphism 
covariant fashion. If the matter field is a gauge field with compact gauge group, we can quantize 
it with exactly the same methods that are used in QGR for the gravitational field. This way, we 
obtain a neat unified description of gravity and the other gauge fields. Also for fermions or scalar 
fields, a representation should be used that is background independent. This rules out the usual 
Fock representation. New representations for fermionic, Higgs and scalar fields in keeping with the 
principles of QGR were proposed in |]25| and we will use them for our purpose. 

The quantization of the Hamiltonian of the coupled system is a rather nontrivial task, due to 
its complicated non-polynomial dependence on the basic variables of the theory. Nevertheless, a 
scheme for the quantization for densities of weight one has been proposed in [|28|, which leads 
to well-defined candidate operators. The resulting operators are quite complicated but perfectly 
well defined and lead to reasonable results in a symmetry reduced context Another very 

encouraging aspect of the scheme is that it works precisely due to the density one character of the 
classical quantities, which is dictated by background independence, and not only despite of it. In 
[m and the present paper we will proceed along the lines given in ||2^, ^ and obtain Hamiltonian 
constraint operators for electromagnetic, scalar, and fermionic fields coupled to gravity. 



As a next step we have to deal with the constraints of the theory: A GauB constraint for Gravity 
and for every matter gauge field, the spatial diffeomorphism constraint of gravity, and, finally and 
most importantly, the Hamilton constraint of the coupled gravity-matter system. 

The implementation of the diffeomorphism constraint has been accomplished in [^]. Still, there 
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is a difficulty related to the spatial diffeomorphism constraint for pure gravity: No spatially diffeo- 
morphism invariant quantum observables (apart from the total volume of the space-like hypersurface 
S, in case it is finite) have been constructed so far. This is due to the fact that such observables 
are given by integrals over S of scalar densities of weight one built from the spatial curvature tensor 
and its spatial covariant derivatives which are highly non-polynomial functions. This problem gets 
alleviated when matter is coupled to the gravitational field. For instance, the matter can serve to 
define submanifolds or regions of S in a diffeomorphism covariant way. Diffeomorphism invariant 
observables can then be obtained by integrating the gravitational fields over these submanifolds or 
regions ||35|, 0. Indeed we will see that this also applies to the Hamiltonian for gravity coupled mat- 
ter: The corresponding operator constructed in the next section will be diffeomorphism invariant. 
This is important for the following reason: Since the diffeomorphisms of S are implemented unitarily 
on Tikin, the expectation value and fluctuations of a diffeomorphism invariant operator do not differ 
from its expectation value and fluctuations in the state that results from projecting the original one 
to the diffeomorphism invariant Hilbert space (via group averaging) ||3^ provided the operator satis- 
fies certain technical conditions (it has to leave cylindrical subspaces of the Hilbert space separately 
invariant). Therefore as long as we work with diffeomorphism invariant operators on Tikin we do 
not have to bother too much about implementing the diffeomorphism constraint. Similar remarks 
concern the GauB constraints, so we will also not be concerned with their implementation in what 
follows. 

We now turn to the implementation of the Hamilton constraint. Even for pure gravity, this is a 
very difficult topic. Though solutions have been found |2^, they are notoriously hard to interpret 
due to the lack of Dirac observables invariant under the motions generated by the Hamiltonian 
constraint (even in the presence of matter) and a thorough understanding of the "problem of time" . 
The problem of finding solutions to the Hamilton constraint for gravity coupled to matter has not 
been treated before although the method of can in principle be applied as well. [| 

Since one of our goals is to explore the semiclassical limit of QGR coupled to matter, the task 
presented to us is even harder: Not only do we have to find some solutions to the Hamilton constraint, 
but we are interested in specific solutions in which the gravitational field is in a state close to some 
given classical geometry. 

As already explained in the introduction, in the light of these difficulties, we propose to proceed 
along slightly different lines. To give an idea what we are aiming at, imagine we ought to compute 
corrections to the interaction of some quantum system (an atom, say) with an electromagnetic field, 
which are due to the quantum nature of the electromagnetic field. Ultimately this is a problem in 
quantum electrodynamics and therefore certainly not solvable in full generality. What can be done? 
For the free Maxwell field, there is a family of states describing configurations of the quantum field 
close to classical ones, the coherent states: Expectation values for field operators yield the classical 
values and the quantum mechanical uncertainties are minimal in a specific sense. Such states could 
be used to model the classical electromagnetic field. Certainly these coherent states are no viable 

^Notice that, since presently the correctness of the classical limit of the operators corresponding to the quantization 
of the geometry and matter Hamiltonian constraints proposed in ^ is not yet confirmed, in order to verify this 
proposal it is well motivated to work with kinematical semiclassical states: This is because one cannot study the 
semiclassical limit of an operator on its kernel. Also, since the spatially diffeomorphism invariant states are not left 
invariant by the Hamiltonian constraint, we cannot even work at the spatially diffeomorphism invariant level. In this 
paper we are, however, not so much interested in testing the Hamiltonian constraint but rather we suppose that some 
correct version of it exists and ask how physical predictions can be extracted without solving the complicated theory 
exactly. 
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states for the full quantum electrodynamics treatment in any sense. They do not know anything 
about the dynamics of the full theory. The key point now is that though being in some sense 
"kinematical" , the coherent states for the Maxwell field are nevertheless a very good starting point 
to compute approximate quantum corrections as testified by the computations in the framework of 
quantum optics [^. Certainly this analogy is not complete in that QED is equipped with a true 
Hamiltonian (rather than just a Hamiltonian constraint) but it shows nevertheless that sometimes 
kinematical states lead to rather good approximations. 

In the present work we will proceed in the same spirit: We will not seek states which are solutions 
to the constraint and approximately correspond to some classical geometry, but rather start by 
considering kinematical semiclassical states. 

Consequently, we treat the Hamiltonians of the matter fields not as pieces of the Hamiltonian 
constraint but rather as observables. In particular, there will be no lapse function, these Hamiltonian 
operators are simply different operators (from the Hamiltonian constraint operator). Although we 
will follow essentially the steps performed in |2^, the fact that we are dealing with different operators 
allows us to change the quantization procedure slightly, for instance, the Hamiltonian operators leave 
the cylindrical subspaces of TCkin separately invariant. 

It is hard to judge the validity of this approach as compared to the desirable full fledged solution 
of the Hamilton constraint. To shed some light on this issue, in an appendix we consider a simple 
quantum mechanical model system. For this system we can show that the expectation values of Dirac 
observables in coherent states states on the kinematical level numerically differ from the results of a 
treatment using dynamical coherent states, the differences are tiny, however, as long as the energy 
of the system is macroscopic. 

Another issue raised by the treatment outlined above is that much depends on the choice of the 
state that is employed to play the role of the semiclassical state. We will defer a discussion of this 
fascinating topic to the companion paper and only make some brief remarks here: 

All candidate semiclassical states proposed so far are graph based states, i.e. cylindrical functions 
in Hkin- Consequently, this is assumed to be the case in the present work. The picture might however 
change substantially if ideas such as the averaging over infinitely many graph based states advocated 
in P could be employed. For some more discussion on this point we refer to [|T], p8| . 



Having adopted the above viewpoint on the Hamilton constraint, we can construct approximate 
n— particle Fock states propagating on fluctuating quantum spacetimes as follows: Denote by 
the gravitational phase space of initial data on the hypersurface S of the differentiable manifold 
M = R X S. Let m G be initial data for some background spacetime. An ordinary n— particle 
state is an excitation of a vacuum state f^matter("^) ^ usual Fock space 'H^^^^^j.{m) which is of a 
completely different type than the background independent Hilbert space T^matter- The construction 
of that vacuum state (and the entire Fock space) makes heavy use of the background metric in 
question, here indicated by the explicit dependence of the state on the point m in the gravitational 
phase space. This dependence slips in because the state f^matterl"^) is usually chosen as the ground 
state of some Hamiltonian operator ffmatter("^) on the background spacetime in question. We now 
see what will heuristically happen when we start switching on the gravitational fluctuations as well: 
The dependence of ^matteri''^) ^ ^'^ become operator valued] In other words, the vacuum 
state function m i— >• ninatter("^) becomes a vacuum operator VL := f2mattcr(^)5 that is, a function of 
the matter degrees of freedom with values in CiTC^^-^) "^^matter where T^matter now necessarily 
a background independent matter Hilbert space (which we have chosen to be the one currently 
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in use in QGR) and £(7-^g™^) denotes the space of linear operators on a background independent 
geometry Hilbert space. We see that the whole concept of an n— particle state becomes a very 
different one in the background independent context! Of course, we do not want a vacuum operator 
but a vacuum state on the full Hilbert space T-^"^'" so that one will apply the vacuum operator to a 
state V'grav('Ti) G ?^grav5 ^^at is, Q{m) = fl'lpm- 

We conclude that a fundamental n— particle state of some matter type corresponding to an ordi- 
nary n— particle state of the same matter type propagating on some background spacetime described 
by the point m in the gravitational phase space will be a complicated linear combination of states 
of the form -ipgravim) ^'matter G T^gJav ® '^matter • How should this statc bc obtained from first prin- 
ciples? We propose the following strategy: Consider the full gravity coupled Hamiltonian operator 
H and construct a annihilation operator from it, which is now an operator on the full Hilbert space 
-^km whose partial classical limit at the point m of the gravitational phase space with respect 
to the gravitational degrees of freedom mirrors the usual Fock space annihilation operator on the 
background spacetime described by m. 

This is what one should do. Now recall that the construction of Fock space annihilation operators 
on a given background involves, for instance, the construction of fractional powers of the Laplacian 
operator on that background metric which is an operator in the one-particle (or first quantized 
matter) Hilbert space. Thus, our fundamental annihilation operator will involve a quantization of 
these Laplacian operators which therefore become an operator on the tensor product of the one 
particle Hilbert space and the gravitational Hilbert space TigJ-av- While we are able to actually 
construct these operators in the present paper, as one can imagine, the formulas that we obtain are 
too complicated in order to do practical computations with present mathematical technology because 
fractional powers of the Laplacian are defined via its spectral resolution which is difficult to find. 

As an approximation to this exact computation we therefore propose to first compute the expec- 
tation value of the Laplacian operator in a gravitational coherent state and then to take its fractional 
powers. Now, precisely because we are using coherent states, this approximation will coincide with 
the exact calculation to zeroth order in h while for higher orders we presently do not know how 
significantly results are changed. The details of these statements will be presented in section |[ 

The last step in the program is then to obtain, in principle testable, predictions from the theory 
obtained so far. For instance, we are interested in states of the form '^/'grav (^) ®'?/'matter ^ '^gmv® '^matter 
and wish then to construct an effective matter Hamiltonian operator as a quadratic form through 
the formula 

< ^matter, -f^m^ttcr('"^) ^matter ^Wj^'^ttcr V'grav("T') ® V'matter, -f^V'grav("^) ® V'matter ^Wljlv^^mattor 

The operator -f^J^tter("^) already contains information about the quantum fluctuations of geometry. 
The quantum fluctuations of matter are certainly much larger than those of geometry in the energy 
range of interest to us, however, as we are interested in Poincare invariance violating effects, which 
are excluded by definition in ordinary QFT on Minkowski space, in order to study those we can 
neglect the quantum effects of matter as a first approximation (that is, we are dealing with free field 
theories except for the coupling to the gravitational field). Therefore, we take the classical limit of 
-^matter l'^) ^nd study the wave like solutions of the matter dynamics it generates. One can also take 
the point of view that this procedure corresponds to the first quantization of matter on a fluctuating 
spacetime. Second quantization will then be studied later on when we discuss n— particle states. 
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Adopting this viewpoint, as soon as a semiclassical state for the gravity sector is chosen, transla- 
tion and rotation symmetry is heavily broken on short scales due to the discreteness of the underlying 
graph. The theory describes fields propagating on random lattices, bearing a remarkable similarity 
to models considered in lattice gauge theory [^, Due to the lack of symmetry on short scales. 



notions such as plane waves and hence dispersion relations can at best be defined in some large scale 
or low energy limit. We will show that the problem of treating these limits is by no means trivial and 
requires careful physical considerations. It is closely related to the condensed matter physics problem 
of computing macroscopic parameters of an amorphous (i.e. locally anisotropic and inhomogeneous) 
solid from the parameters of its microscopic structure. 

To get a feeling for the problem, we will sketch a one dimensional model system for which we are 
able to find exact solutions. We will then turn to general fields on random lattices and describe a 
procedure to obtain dispersion relations valid in the long wavelength regime. 

This concludes our explanatory exposition. We will now proceed to the details. 



3 Review of Quantum Kinematics of QGR 

In QGR the manifold underlying spacetime is taken to be diffeomorphic with M = M x E where E 
represents a 3d manifold of arbitrary topology. We will now summarize the essential aspects of the 
kinematical framework of pSf for matter fields coupled to quantum gravity. We also introduce the 



Hamiltonians that we will be deriving dispersion relations for. 
Gravity and Gauge Theory Sector 

The canonical pair consists of a G connection Al^{x) for a compact gauge group G and a Lie{G) 
valued densitized vector field Ef{x) on E. Here we can treat all four interactions on equal footing. 
For the gravitational sector we have G = SU (2) and the relation of the canonical pair to the classical 
ADM variables qab, Kab is 

dei{q)q'^' = cEtE^^, = - ^^K^bE^", 

Vdet(g) 

where F is the spin connection corresponding to the triad E, and t is the Barbero-Immirzi parameter 



which can in principle take any nonzero value in C [0, We will choose l = 1 in what follows. 
As for units, we choose [A] = meter"^. As a consequence, E will be dimensionless for gravity and 
has dimension cm~^ for Yang-Mills theories. 

In the following we will have frequent opportunity to use the notion of graphs embedded in E: 

Definition 3.1. By an edge einJ^we shall mean an equivalence class of analytic maps [0, 1] — > E, 
where two such maps are equivalent if they differ by an orientation preserving reparametrization. 
A graph in E is defined to be a set of edges such that two distinct ones intersect at most in their 
endpoints. 

There is some notation in connection to graphs that we will use frequently: 
The endpoints of an edge e will be called vertices and denoted by 6(e) (the beginning point of e), 
/(e) (the final point of e). 
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The set of edges of a graph 7 will be denoted by E{'~f), the set of vertices of its edges by V{'y). 
Given a graph 7, we will denote the edges of 7 having v as vertex by £"(7, v) or E{v) if it is clear 
which graph we are referring to. 

Given a graph 7, a vertex v e ^^(7) and an edge e e E{v) we define 

if 6(e) = V 
a{v,e) = {-I if /(e) = v 

if e is not adjacent to v 

Thus e'^'^'"'^^ is always outgoing with respect to v. 

Being a one-form, A can be integrated naturally (that is, without recurse to background structure) 
along piecewise analytic curves e in S, to form holonomies 



he[A] = Pexp 



A 



e G. 



It is convenient to consider a class of functionals of the connection A a bit more general: 

Definition 3.2. A functional f[A] of the connection is called cylindrical with respect to a piecewise 
analytical graph 7 if there is a function 

f : G-I^WI — , C, 

such that 

f[A] = f{h,,[A],h,,[Al ...), €,,62,... e E{j). (3.1) 

The density weight of E on the other hand is such that, using an additional real internal vector 
field /* it can be naturally integrated over surfaces S to form a quantity 



Esj = [ r {*E), 
Js 



analogous to the electric fiux through S. 

In the connection representation of diffeomorphism invariant gauge field theory, quantization of 
the Poisson algebra generated by the classical functions Cyl and the vector space space of electric 
fluxes £ is achieved on the Ashtekar-Lewandowski Hilbert space 

It is based on the compact Hausdorff space A of generalized connections which is a suitable enlarge- 
ment of the space of smooth connections A and the uniform measure //q- 
The classical Yang-Mills Hamiltonian (coupled to gravity) reads 



Hym = / d'x^^[E'^,E'j + B'iB'j]5''^. (3.2) 

^VyM JS ydet^^j 

Here Qym is the Yang-Mills coupling constant, E'} is the Yang- Mills electric field, i?j = e^'^^F^^ the 
magnetic field associated with the Yang-Mills curvature Fj^. 
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We now quantize this operator along the hnes of , actually only for Maxwell Theory since in 
this paper we are interested only in free theories when taking the metric as a background field. Here 
we take advantage of the fact that Hym is an operator of its own although, of course, the integrand of 
( |3.2| ) is a piece of the classical Hamiltonian constraint of geometry and matter. Accordingly we may 
exploit the quantization ambiguity concerning the loop attachment in as follows: We define the 
operator Hym consistently on the combined spin-network basis of matter and geometry introduced 
p5[] and use the following notion. 



m 



Definition 3.3. Let a graph 7, a vertex v G V{'-^) and two different edges e, e' G i?(7) incident at v 
he given. By a minimal loop based at v we mean a loop I3{'j,v,e,e') in 7 which 

• starts at v along e and ends at v along e' , 

• does not self-overlap, 

• the number of edges used by [3 except e, e' cannot be reduced without breaking the loop into 
pieces. 

Notice that given 7, f , e, e' a minimal loop does not need to be unique! Denote by S'(7, f , e, e') the 
set of minimal loops corresponding to the data indicated and by L(7,f,e,e') their number. Notice 
also that the notion of a minimal loop does not make any reference to a background metric, it is an 



object that belongs to the field of algebraic graph theory |4^, ^ . 

Here we see the first difference as compared to the loop choice in ||2^: A minimal loop is always 
contained in the graph that we are dealing with. The second difference that we will introduce in 
contrast to is that there we used functions of holonomies of the type Hp — H^^ in order to express 



the Yang-Mills magnetic field in terms of holonomies. However, as correctly pointed out in the 
regularization ambiguity allows more general functions, the only criterion is that the final operator 
is gauge invariant and that the function should vanish at trivial holonomy. Our preliminary proposal 
for the Maxwell Hamiltonian operator, projected to spin-network states over graphs 7 is then 



1 r 1 

l,e2-Le 



X {-Y,Y,, + — — [ V — = V 

P(7, f , e)P(7, f , e'jo; ^-^ L(7,f, 61,62) i 



xl E r, - M E — ]}. 

^-^ L7,t>,6i,6, ^-^ I 

where for a vertex w , a real positive number r and an edge e starting at it we have defined the basic 
operator 

Qi{v,r) := ^^iT{rMK\{V,r]). (3-4) 

with Vu the volume operator over an arbitrarily small open region containing v. We are 

using a basis of su{2) with ii{TjTk) = —'26jk- Here a = hQM is the Feinstruktur constant and mp, Ip 



are the Planck mass and length respectively. We have adapted the coefficients of [27| to the case of 
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G = U{1) and we have distinguished the Maxwell holonomy H from the gravitational holonomy h. 
Ye is the right invariant vector field on U{1) with respect to the degree of freedom Hf.. The notation 
e _L e' means that e 7^ e', e fl e' = f 7^ and that the tangents of e, e' at v are linearly independent. 
The number P(7, v, e) is the number of pairs of edges ei, 62 with t> G ei fl 62; ei, 62 J- e and A^(7, v) 
denotes the valence of the vertex v. For each vertex v the edges incident at v are supposed to be 
outgoing from v, otherwise replace e by e'^^"'^-' everywhere in ( p.3|) . The branch of the logarithm 
involved in ( p.3|) is defined by ln(l) = 0. The logarithm is convenient in order to define photon states 
later on but any other choice will do as well, just giving rise to more quantum corrections. 

The manifestly gauge invariant and spatially diffeomorphism invariant Hamiltonian ( |3.3| ) is pre- 
liminary because we may want to order it differently later on. Notice that it does not have the correct 
classical limit on an arbitrary graph, the graph has to be sufficiently fine in order to reach it! We 
will show that it defines a positive definite, essentially self-adjoint operator on Tikm- 

Lemma 3.1. For any positive real r the operator iQi{v,r) on Tio defined by ( ^.4 ) is essentially 
self-adjoint with core given by the core ofV^. 

Proof. Since {hejAB is a bounded operator it suffices to show that iQi{v,r) is symmetric with dense 
domain the core of Vy. 

Using that [{he)AB]^ = {K^)ba and {tj)ab = -{tj)ba we find 

At[i Qiiv, r)]t = {(^ [{ih:')cAy , V:] {{K)bc)^ = -i{r,)BA [((Mac) , KT] {{K')cb) 
-i Tr I Tj 



Tr(rA [h~\V:])=Ar[iQi{v,r)] 



because Tr(rj)K' = 0. □ 
Scalar and Higgs Fields 

We will consider only Lie(G) valued Higgs fields 0* with canonically conjugate momentum tTj. In 
particular, a neutral scalar field is Lie(t/(1)) valued and transforms in the trivial adjoint represen- 
tation. We will take 0i to be dimensionless, then vTj oc tt* has dimension cm~^. 

The background independent Hilbert space of |^ is based on the quantities 



U{x) := e'^'(^)"' and tirj = / d^xf'ni, (3.5) 



R 

where U{x) is referred to as point holonomy, tj is a basis of Lie(G) and the second quantity is 
diffeomorphism covariant since Hi is a scalar density. One can then quantize the Poisson algebra 
generated by these objects on a Hilbert space L2(U, djjLu) where W is a distributional extension of the 
space lA of smooth point holonomies and d^u is an associated uniform measure. This Hilbert space is 
very similar in spirit to the one for gauge theories displayed in (P). A dense subspace of functions in 
this Hilbert space consists of the cylindrical functions. Here a function is cylindrical over a graph 7, 
if it depends only on the point holonomies f/(f ), v G ^(7)- That the point holonomies are restricted 
to the vertices of a graph is dictated by gauge invariance (for neutral scalar fields there is clearly no 
such argument but given a function depending on a finite number of point holonomies we can always 
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trivially extend it to depend trivially on gravitational holonomies over a graph with the arguments 
of the point holonomies as vertices). 

In this paper we are only interested in neutral Klein-Gordon fields without interaction potential. 
The unitary operator U{x) acts by multiplication while the momentum operator is densely defined 

by 

veV{'y)r\R 

where Y^j denotes the right invariant vector field on U{1) with respect to the degree of freedom U{v) 
and Qkg is the Klein Gordon coupling constant which is such that HQkg has dimension cm^. 
The classical Klein Gordon Hamiltonian coupled to geometry is given by 

Hkg = 7^ / d'x[ Z— + v/dit(^[g"V,a0,. + K'<P% (3.6) 
^(^KG Jt. A/det(g) 

where is the Compton wave length of the Klein Gordon field. To quantize ( |3.6| ) we again copy 
the procedure of and define it on combined matter - geometry spin-network states with the 
following modifications: 1) no new Higgs vertices on edges of the graph are introduced and 2) we 
replace the function \U{x) — that substitutes (j){x) by something else in accordance to what we 
have said for gauge fields already. We then propose the preliminary version of the Klein Gordon 
Hamiltonian operator, projected to spin-network states over graphs 7 by 

^KGrf = — ^^^P 1^ ^ 
^ DgV(7) 

' ?jgeine2ne3;ei±e2-Le3 

V e'^^e [MW(e/))-ln([/(b(e/))] ^. 3 3 ^ 

W ^ V e^M^e NW(eL))-ln([/(6(e.))] .^ 3^„ 3 

^ J5|)L„„, 5^ MtM Mt^lK,. ,3.7) 

fSv (7) 

where T(7, v) is the number of triples of edges incident at v with linearly independent tangents there 
and 6(e) and /(e) respectively denote starting point and end point of an edge. The operator ( p. 7] ) is 
again manifestly gauge and diffeomorphism invariant. Notice that hQxG/^'p is dimensionless while 
has dimension cm^ so that all terms have mass dimension. We already have ordered terms in (|3.71 ) 
in a manifestly positive way and the branch of the logarithm used corresponds to the fundamental 



X 



X 
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domain of (D again, i.e. ln(z/|2;|) G {—n,n]. 



Fermion Fields 



Since the canonical, non-perturbative quantization of the Einstein-Dirac theory in four spacetime 
dimensions using real valued connections is maybe less familiar to the reader we review the essential 



aspects from ||2^ in slightly more detail. The classical Hamiltonian reads (we neglect coupling to the 



Maxwell field in this paper since we want to isolate effects of quantum gravity on the propagation of 



free fields, see |2^, for the coupling to non-gravitational forces) 



H=hj (fx{^[VaJ, + i{i^'^aPa^ - faPaTi ' c.c.) - Ki{^^ ^ - fri)\ + 2iroVdet(g)(^^77 - f^)} 

(3.8) 

where is the rest frame wave number, cxj are the Pauli matrices, Jj := ip'^ajip + ffajT] is the 
fermion current, ip = and 7] = {rjA') respectively denote the left - and right handed components 
of the Dirac spinor = {tp, rj)'^, ip denotes the involution on GraBmann variables and the complex 
conjugation c.c. is meant in this sense. The spinors ip, rj transform as scalars under diffeomorphisms 
and as left and right handed spinors under SL{2,(C). In particular, Vaip = da'p + ^A{Tj'p,T)a'r] = 
daV + ^^i'^jV where tj = —icrj. Our convention for the Minkowski space Dirac matrices is 7° = 
—ia2 ® I2, = 0"! ® <Jj appropriate for signature (— , +, +, +). The dimension of our spinor fields 
is cm~^/^ so that ( p.8|) has indeed dimension of energy. Notice the explicit appearance of the field 
I{J = - P . 

a a a 

A peculiarity of spinor fields is that they are their own canonical conjugates. Consider the half- 
densities 



e := ^det(g)^, p := ^det{q)r], (3.9) 
then the canonical anti-brackets are given by 

{U{x),pB{y)}+ = , {pA{x),pB{y)}+ = , (3.10) 



while all other anti-brackets vanish. That ( |3.9D mixes gravitational and spinor degrees of freedom is 
absolutely crucial: without this peculiar mixture it would not be A{ that is canonically conjugate 
to Ej/k but rather A{ + iipei{'p'^'p + fj^rj) which is now complex valued and this would destroy the 
Ashtekar Lewandowski Hilbert space Tio since connections would become complex valued. 

Clearly, we obtain the Einstein- Weyl Hamiltonian by setting either ^ or p to zero. Likewise we 
can treat the case of several fermion species by adding appropriate similar terms to ( |3.8| ). In what 
follows we just stick with (|3.8| ), the reader may introduce the appropriate changes for the case by 
hand himself. 

In order to quantize (|3.8|) we want to write ( p.8|) into a more suggestive form. To that end, notice 
that the Gauss constraint in the presence of fermions reads 

1 E^ h ^ ^ 
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where k is the gravitational constant and where T>a acts on tensorial indices by the Christoffel 
connection associated with g^b and on SU{2) indices by the connection A^. Thus we can solve the 
Gauss constraint for the fermion current so that after an integration by parts we have the identity 

n f d'a:%VaJ,] = - f d'x^^^£!L (3.12) 



2 ' « is v/dit(^ 

modulo the Gauss constraint which now depends only on the gravitational degrees of freedom. 
Clearly, in flat space ( |3.12 ) vanishes. 
Next it is easy to see that 



El 
Vdet(g)' 



+zE]{2P^ajVatp - faPaV - c.c.) = % ^-l—- {l^aPai - P^f^PaP - c.c), (3.13) 



where Vai ■= dai + v4^rj,^/2 ignores the density weight of (and similar for p) since the appropriate 
correction term is cancelled through a similar term in the c.c. piece. 

Formulas (|3.12| ) and ( p.l3|) imply that (|3.8|) can be rewritten in terms of as 



1 f ,, {VaE. 



a\2 



H = - I d'x '' " (3.14) 

+ n / Sx{-^==\-Vi{l^ap,i - p^ap^p - c.c.) - Kii^i - p^ p)] + iK^ii" p - p^i)}. 
Jt, 2-^det(g) 

(3.15) 

Now recall from p5| in that for reasons of diffeomorphism covariance it turned out to be crucial 
to work instead of with the half densities ^, p with the scalars 

e^ix) := f d'y^6{x,y)U{y) , 0'^{x) := [ d^y^6{x,y)pA{y), (3.16) 



which still transforms covariantly under gauge transformations 6{x) — > g{x)6{x) since the distribution 
^y 6{x, y) has support aX x = y. If we require the fields 9 to be ordinary GraBmann fields, then 
formula (|3.16|) implies that the spinor half-densities ^, p are distributional GraBmann fields. This 
distributional character is due to the factor '^det(g) which in quantum theory becomes an operator 
valued distribution proportional to \j5{x, y) (recall that there is no such thing as classical fermion 
fields). The inversion of ( p.l6|) is given by 

U^) := ^5{^eA{y) , Pa{x) := Yl VSi^O'M, (3.17) 
2/es yes 

due to the identity a/ S{x, y)S{x, z) = 6{x,y)6y^z where 6x,y denotes the Kronecker symbol (equal to 
one when x = y and zero otherwise). 

Let fe{x,y) = fe{y,x) = fe{x — y) he Si one parameter family of smooth, nowhere negative 
functions of rapid decrease such that also \/Jj^x~y) is smooth and such that lim^^o fe{x, y) = 5(x, y). 
An example would be f,{x,y) = Ua [e'^^^-^^^'/^^^Vv^]. Then 

{dj){x) := hm / d'y{d,a./fj^)^{y) 
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-lim / d'y{dya./fj^))^{y) 



s 



= lim / d'y./fj^){dai){y) 

= jjhVK^{dai){y) (3.18) 

where in the integration by parts no boundary term was picked up since /g is of rapid decrease. It 
follows that 

{VJ){x) = [ d?y^/5{^){V,i){y) ^ = E (I^a^) (y) (3.19) 

y 

for classical (smooth) A^. 

The Fermion Hilbert space now is constructed by means of Berezin integral techniques where our 
basic degrees of freedom are the 9a{x), 9'^{x) and their involutions. We start with only one GraBmann 
degree of freedom and denote by S superspace with anticommuting GraBmann coordinates 9, 9, that 
is, 9^ = P = 0, 99 = -99. 

A "holomorphic" function depends only on 9 and not on 9 and is of the general form 

f(9) ^a + b9 (3.20) 

with arbitrary complex valued coefficients a, b while a generic function on S is of the general form 

F(9, 9)^a + b9 + c9 + d99 (3.21) 

with arbitrary complex valued coefficients a, b, c, d. The integral of F over S with respect to the 
"measure" d9d9 is given by 

J d9d9F{9, 9) = d. (3.22) 

A quantization of the canonical anti-brackets 

{9, 9}+ = {I ~9}+ = 0, {1 9}+ - {9, 9}+ = ^ (3.23) 
and of the reality conditions 

0^e, 9 = 9 (3.24) 

can be given on the space L2{S,diJ,F) of "square-integrable" holomorphic functions with respect to 
the "probability measure" 

duF := e^^ d9 d9 = [l + 99] d9 d9 (3.25) 

which is positive definite: 

< /, /' >:= J dfiF{9, 9)7(9) f (9) = da' + W for f{9) = a + b9, f'{9) = a' + b'9. (3.26) 
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We just need to define the operators 0,9 hj 

{Ofm := Of (6) = aO, {Ofm ■= ^/W = b (3.27) 
(derivative from left) and verify immediately that the canonical anticommutation relations 

[§j]^ = 2e^ = [le]+ = 2f = 0, [e,e]+ = [e,e]+ = ee + fe = i (3.28) 

as well as the adjointness relations 

= e, f = e (3.29) 

hold with respect to the measure dfip- 

This covers the quantum mechanical case. Let us now come to the case at hand. Recall that we 
had the following anti-brackets for our spinor degrees of freedom 

{Ux),Uy)}+ = {U^),Uy)}+ = 0, {U^),Uy)U = {Uy),U^)U = (3.30) 

and similar for p. Inserting the transformation ( p.l6| ) we see that ( |3.30| ) is equivalent with 

{^A(x),^B(y)}+ = {eAix),eBiy)}-, = 0, {eAix),eBiy)}+ = {eBiy),0Aix)}+ = (3.31) 

so the 6 distribution is simply replaced by the Kronecker symbol. This suggests to define the Fermion 
Hilbert space as the continuous infinite tensor product ||5^ 

L2{S,dfiF) (3.32) 
where 9a{x),9a{x) = Oa^xY are densely defined on Co vectors by 

9A{x)^f := [<^x^y,Bfy,B] ® [fx~A ® (^/x,a)], 

OAixY^f := [/,,_a®(^Vx.,a)]. (3.33) 

This Hilbert space is unnecessarily large for the following reason: Due to gauge invariance the spinor 
fields are confined to the vertices of an at most countably infinite graph. In particular, if we are 
dealing with finite graphs only, then the subspace Hd of the Hilbert space (|3.32|) , defined as the 
inductive limit of the cylindrical spaces 

H^^D '■= ®»;GV{7),A=±l/2-^2('5, (i/iir) (3.34) 

via the isometric monomorphisms U^y for 7 C 7', densely by 

U^y : Hy^D ^ ^7',d; ®/ := ^v€Vi-y),Afv,A ^ [®v£V{-y),Afv,A] ['^x€V{'y')-V{^),B'^] (3.35) 

is completely sufficient for our purposes in this paper (as long as cr is compact, otherwise we can 
use the techniques from |Q). Equation ( |3.35| ) displays TYd as the strong equivalence class Hilbert 
subspace of T-C% formed by the Cq vector 1 := ^x,a^- 
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We now turn to the quantization of ( p.l4 ). Actually we will not consider the terms which are 



proportional to K^, (VaEj) because they vanish in flat space (with which we are mainly concerned 



in this paper). Of course, quantum corrections will give non- vanishing corrections but since we are 
doing only exploratory calculations in this paper, let us just not discuss those terms. Then the 
methods of lead to the the following quantum operator restricted to matter - geometry spin 
network functions over a graph 7 



' i;Geine2ne3;ei±e2-Le3 
' D'6eine2ne3;ei±e2ie3 

-ihKo E SABS.A^>')0\iv) - enivy'liv)]. (3.36) 

v,v'GV{-y) 

It is not difficult to see that this operator is self-adjoint. Again, as compared to we have chosen 
a different ordering and there are no new fermion vertices created on cylindrical functions over 7. 

Remark 

We have defined the operators Hm, Hkg, He in the combined spin-network basis of matter 



and geometry defined in ||25[. Such a spin network function Ts{A, Am,4',0,6') carries a label s = 
(7, j, n, rh, B, B', I) consisting of a graph 7, a coloring of its edges e with Einstein non-zero spins je € j 
and non- zero Maxwell charges G n as well as a coloring of its vertices v by non-zero scalar charges 
G m, non-zero left-handed fermion helicities B^ G B, non-zero right-handed fermion helicities 
B'^ G B' and intertwiners /t, G / which make the state gauge invariant under the action of the gauge 
group SU{2) X U{1). This defines densely a continuum operator by HTg := H^(^s)Ts on the Hilbert 
space H = He ® Hm ® Hkg ® 'Hd {E stands for the Einstein sector) where H = Hm + Hkg + Hd. 
However, the operators are not the cylindrical projections of H since they are not cylindrically 
consistent, i.e. H^Cy\^, 7^ H^iCyl^, for 7' C 7. Rather, in order to evaluate the operator on cylindri- 
cal functions, one has to decompose them in terms of spin-network functions. Nevertheless one can 
construct from the family {H^) a cylindrically consistent family {H"') as follows: 

Notice that each has the following structure 

-^7 = E -^'T''"' (3.37) 
t)ev(7) 

where is a local operator, that is, it depends only on the finite subset E^i^'j) C -£^(7) of edges 
of 7 incident at v. For e G -^(7) denote by Pg the projection operator on the closed linear span 
of cylindrical functions in Cyl^ which depend through non-zero spin on the edge e. For subsets 
EcE' C E{-f) let 

W':=[nA.e][ n (1-A,e')]- (3.38) 

eS-B e'eE'-E 
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Consider now the operator 



■- ^ P-y,E,E^i'r)H^~[E^('y)-E],vP-y,E,E^i-,) (3.39) 
Sc£i,(7) 

where the sum is over the power set of E^^'j) (set of all subsets) and with it the, still self-adjoint, 
cylindrically consistent family of operators 

H"i = J2 HI (3.40) 

Notice that H^_[E^(^yyE],v = whenever \E\ < 3. 

It is important to notice that both families {H"^) and (-^7) give rise to the same continuum 
operator as long as there are no gravitational holonomy operators outside of commutators involved 
(as is the case for the bosonic pieces). It is just sometimes more convenient to have a consistent 
operator family if one does not want to decompose a cylindrical function into spin-network functions. 
In fact, our semiclassical states over 7 are not spin-network states over 7 but they are cylindrical 
functions, i.e. linear combinations of spin-network states where also all smaller graphs 7' C 7 appear. 

The careful reader will rightfully ask whether the expectation values of the operators defined in 
terms of ( |3.37D and ( ^.401 ) respectively agree on semiclassical states over 7 (they do exactly if no 
gravitational holonomy operator is involved). This is important since the operator P'y,E,Ey(-y) did not 
come out of the derivation in |^ (for operators involving the gravitational holonomy) and thus could 



spoil the classical limit. Fortunately, the answer to the question is affirmative due to two reasons: 

1) The expectation value of P^^e turns out to be of the form 1 — e"'^/*'^ where c, P are positive numbers 
of order unity for non-degenerate metrics and t is a tiny number related to h. Then the expectation 
value of 1 — P^ e is of order 0{t°°). 

2) We will use only graphs for semiclassical calculations such that the valence of the vertices 
is bounded from above. Thus the sum ( p.39|) involves only a small number of terms and (1 — 



g-c/i^)|E.(7)l = l+0(t°°). 

Thus the expectation value of with respect to cylindrical functions over 7 agrees with that of H'^ 
to any finite order in t and for semiclassical calculations we can practically treat the family (i^^) as 
if it was cylindrically consistently defined. We will assume that to be the case in what follows. 



4 Matter n— Particle States on Graphs 

The aim of this section is to sketch how one would in principle construct exact n— particle states 
propagating on fiuctuating quantum geometries as well approximations of those by using quantum 
geometry expectation values in gravitational coherent states. The computation of those expectation 
values is sketched in the next section, more details can be found in our companion paper Since 
in these two papers we are only interested in qualitative features we do not want to spend too much 
technical effort and therefore make our life simple by replacing SU{2) by U{1)^, see [|10|, |lT| for how 
non-Abelian gauge groups blow up the computational effort by an order of magnitude. We leave the 
exact computation for future investigations after the conceptual issues discussed in this work have 
been settled.0 

^Actually, SU{2) is replaced by U{1)^ in the GNewton limit if one rescales the gravitational connection A by 
A/ G Newton (lonii-Wigner contraction), but GNewton — *■ also implies Ip ^ Q and this is precisely the regime we are 
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4.1 Specialization to Cubic Random Graphs 

A second simplification that we will make is to consider for the remainder of this and the companion 
paper only graphs of cubic topology. Graphs of different topology can be treated in principle by the 
same methods that we develop below but for analytical computations graphs of different topology 
present an extremely hard book keeping problem which is presumably only controllable on a com- 
puter. But apart from these more practical considerations we can also give some physical motivation: 
A) 

As is well-known, there exist an infinite number of discretizations of a classical continuum action or 
Hamiltonian with the correct continuum limit and some of them reflect the continuum properties of 
the action or Hamiltonian better than others. In that respect it is relevant to mention the existence 
of so-called perfect actions |^T| which arise as fixed points of the renormalization group flow for Eu- 
clidean field theories. These are perfect in the sense that although one works at the discretized level, 
the expectation values are, for instance, Euclidean invariant] These techniques have been applied 
also to differential operators and there exist, for instance, Euclidean invariant Laplace operators on 
arbitrarily coarse cubic lattices |^ |5^ despite the fact that a cubic lattice seems to introduce an 



unwanted direction dependence! The quantization of our Hamiltonian operators could exploit that 
freedom in order to improve semiclassical properties, the choice that we have made is merely a first 
natural guess. 
B) 

Secondly, the graph does not need to be "regular" but rather could be an oriented random cubic 
graph adapted to the three metric qat to be approximated. In fact, we will discuss this possibility in 
detail in our companion paper when we discuss (light) propagation on random cubic graphs. Such a 
graph could be obtained by a suitable random process. Let us sketch a procedure for two dimensional 
Euclidean space: We start with a sequence of randomly chosen vectors Vi,V2, ■ ■ ■ subject to the con- 
dition that the angle between vectors adjacent in the sequence lies in the range [—it/2, 7r/2]. Then a 
first sequence of vertices of the graph could be obtained as 0, -iTi, -|- • • • and a first sequence of 
edges by connecting the vertices by straight lines. Now choose another sequence of vectors subject to 
the same condition on the angles, as well as an additional vector p. Again we obtain a set of vertices 
p,p + Vi,p + Vi+V2, ■ ■ ■ and corresponding edges. These are to be discarded if an edge intersects one 
of the edges obtained before. Otherwise the vertices and edges are added to the graph. This can be 
iterated until a convenient number of vertices has been obtained. Then the rest of the edges neces- 
sary to turn the graph into a cubic one can be obtained by connecting the vertices "vertically" by 
straight lines. This description is certainly sketchy, but it could easily be made precise by supplying 
the details of the random distributions for the choice of the points, directions etc. Similarly, it can 
be generalized to three dimensions and curvilinear edges. 

Although it would be extremely cumbersome to obtain analytical results about the resulting 
random graphs, procedures like the one sketched above can be implemented on a computer in a 
straightforward way and any desired information can then be obtained numerically. An artistic 
impression of a small part of a random graph of cubic topology is given in figure |l|. Note that it does 
not favor a direction on a large scale, though it is certainly not rotationally symmetric in a strict 
sense. 

A random cubic graph is certainly still diffeomorphic to a regular cubic graph in with its 
natural Cartesian orientation. Its vertices v can be thought of as points in Z^, its edges can be 



interested in. However, ultimately we must do the SU (2) computation. 
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Figure 1: Part of a (two dimensional) random graph of cubic topology. 



labelled as ej{v), I = 1,2,3 with b{ei{v)) = v, f{ei{v)) = v + I where v + I denotes the next 
neighbor vertex of v along the /—direction. Given a 3-metric Qab to be approximated, let e be the 
average length of the e/(w) with respect to it. Notice that the angles and shapes of the edges are 
completely random! Depending on the random process that has generated the cubic graph, within 
each close to flat coordinate patch there will also be an isotropy scale 6 = Ne at which the the graph 
looks homogeneous and isotropic. The meaning of the scales e, 6 with respect to physical processes 
such as light propagation which introduces a length scale of its own, namely the wave length A, is 
subject to a discussion in our companion paper where we will see how these scales fit together with 
the quantum gravity scale ip and what their relative sizes should be. For the purposes of this section 
we just need that the graph 7 in question has cubic topology. 

We can now specialize our formulas for the family of matter Hamiltonians (H-y) to cubic graphs 
7. Given v G ^(7) let e/(f ) := e/(f ), ej{v) := {ej{v — so that b{ef{v)) = v. Then it is easy 
to check that the local volume operator becomes 




K = 4V ¥'\-^ ^-^] V-^ — 2 — ^^-^^ 

where denotes the right invariant vector field on SU{2) with respect to the degree of freedom h^. 
Notice that we used the coefficient 1/(8 ■ 3!) that was derived in ES, As we will see, with respect 



to coherent states of the type constructed in |1^, |Tl|], only cubic graphs will assign correct expectation 
values to the volume operator if the coefficient 1/(8 ■ 3!) is usedl 

For a graph of cubic topology we have P(7,f,e) = 4, T{'j,v) = 8 since each vertex is 6-valent. 
Also, there are 12 minimal loops based at v along the edges e'^^ (v) , Cj^ (v) , I ^ J each of which 
is unique, that is L(7,t>,e,e') = 1, due to our simple lattice topology. Given v,I,a there are 4 
minimal loops 0-2 (^) ^^^^S the edges ej^{v),e'^{v) with ejj^x = 1 whose orientation we choose 
to be such that the tangents of e'^{v), ej^{v), e^(t') at v in this order form a 3 x 3 matrix of positive 
determinant. With the notation h'^{v) := he'^(v), cr = ±1 and Q-^j^{v,r) := <5e<T(^)('y, t) our matter 
Hamiltonian constraint operators ( |3.3| ), ( p.7| ) and ( p.36| ) become, respectively: 

Maxwell Field: 

^ t)eV(7) /, J,o-,cr' 



20 



IE 



HHpJ^, ,,(.)) 



-]} (4.2) 



a"l,o-2 



where Yj„{v) 
Remark: 

A remark is in order concerning the logarithms that appear in ( [4 .21) . Recall that in the previous 
section we defined for a loop f3 the number ln{Hf^) by the main branch of the logarithm, specifically 
ln{Hj3)/i G [— TT, vr). Now for a minimal loop f3 and sufficiently fine 7 it is indeed true, for a classical 
connection, that \n{Hp) = ln(ifej + .. +ln(ifeiv) if = ei o .. o cat where all appearing logarithms are 
with respect to the same branch. Notice that the right hand side is gauge invariant only for small 
gauge transformations. In the case of a cubic graph we will show below that this relation can be 
written, e.g. ln(iJ^/^ ^(^j) = e^^^^^J \ii{Hejf{v) where {d^f){v) = f{v+I) — f{y). Using the transversal 
projector P-*- of ( [4.16[ ) we can further write this as ln(if^/^ ^(^j) = e^J^ . djP^ ■ ln(ifeK{«)- Finally, 
again for sufficiently fine 7 we may write this as ln(if^/^ ^(^)) = e^'^^ dj[{P-^ ■hi{He^(^y)){mod27i)]. The 
quantity in the square bracket is now manifestly gauge invariant even after extending to distributional 
connections. It is this definition that we will be using for ln(if^i^ ^(„)) in what follows, it simply 
corresponds to a different choice for the quantization ambiguity (for cubic graphs). 

When we compute commutators of such logarithms with electric flux operators on sufficiently 
fine graphs we must regularize the commutator by first restricting to classical connections (so that 
the value of the logarithm lies in (— tt, tt) and then extending the result to A. This being understood 
we get [Y/o-(f ), ln(if/3)] = Hj^'^[Yi(j{v), H13] without picking up discontinuities. Similar remarks hold 
for the scalar field. 



Klein-Gordon Field: 
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where for a function F : ^^(7) hh. (D we have defined the edge derivative [deF]{v) := F{f{e)) — F{v) 
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if V = 6(e). Specialized to the cubic graph we write {dajF){y) := ((9ej(t,)F)(t>) and {djF){v) : = 
a{do-jF){v) is the forward (backward) edge derivative at w if ex = 1 (cr = —1). 
Dirac Field: 

0"liO'2,0"3 

-ihKo Yl SABS^AWW^iv) - eBiv')9'liv)]. (4.4) 

i>,f'Gy(7) 

The operator (|4.4| ) is not manifestly positive definite on a non-fiat background. Fortunately, in con- 
trast to the bosonic Hamiltonians, positivity is not required in order to arrive at suitable annihilation 
operators since it is already normally ordered (subject to the usual particle - antiparticle reinterpre- 
tation upon passage to second quantization). In order to obtain positivity of this Hamiltonian we 
have to invoke a positive and negative energy decomposition of the 1-particle Hilbert space as done 
in QFT on CST. This step will not be preformed here as it goes beyond the exploratory purposes of 
this paper. 

4.2 One-Particle Hilbert Spaces on a Graph 

The operators in (|4.2| ), (|4.3|) and ( |4.4| ) define operators on the Hilbert space Hkm = '^fm ® '^tin ® 
® '^kin with spin-network projections of the form 

H,= Y Mi{v)\{v)G^,^i),^,,,v)Mv{v'), (4.5) 

v,l,v',l' 

where f , v' G ^(7) and /, /' are elements of a discrete label set C of labels like the labels j, /, a, A, fi 
where /x = 1,2 and 6^ = 6,6'^ = 9'. Mi{v) is a linear matter operator for each pair (f , /) while 
G(v,i),(v',i') is a geometry operator for each pair of pairs (f , /), (f', I'). Our aim is to reorder in such 
a way that it acquires the usual form in terms of creation and annihilation operators. The Fermion 
Hamiltonian is already in this desired form, because 

already satisfies the canonical anticommutation relations but the bosonic terms do not. In order to 
do this we need to introduce the one particle Hilbert spaces on the graphs 7. These are defined as 
spaces of complex valued functions F : V{'j) x £ ^ (D; {v,l) 1— > Fi{v) which are square summable, 
that is, their norm with respect to the inner product 

<F,F'>^x:=5^:F)Mi^/(^) (4.6) 

v,l 

converges. Thus, H!y = ^2(^^(7) ® -C) is equipped with a counting measure. Next to this we also 
introduce the Hilbert subspaces .U^^ .U^^ .U^ of square integrable cylindrical functions over 7 
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of the respective field type. Having done this, we can consider the gravitational operator G 
CiTi^) (£(.) denotes the space of hnear operators over (.)) as an operator G G CiTi^ ® Til^) densely 
defined by 

[GTf®F]{A,{vJ)) := Y.^G^^M<v)Tf]{A) F,{v'), (4.7) 

{v',l') 

where Tf is a gravitational spin network state over 7 and A ^ Ae is a. gravitational generalized 
connection. That the right hand side of ( |4.7| ) is indeed again an L2 function will be shown below. 

Likewise, we may consider the operators Mi{v) E C{'H^^^^'^^) as operators M : 7-^™'^**'='' H^^^^'^^x 
Til^ densely defined by 

[MTr"^'-](A,(t;,/)) := [Mi{v)Tf]{A) 

where Tf^^^^^"^ is a matter spin network state. 

But even better than that, for the bosonic pieces of ( [4.6| ) we will be able to show that the operator 
G is positive definitel By inspection then, the whole (bosonic piece of the) operator Hy is a positive 
operator. Moreover, we will be able to take square roots of this operator, defined in terms of its 
spectral resolution on the Hilbert space ® H}^. These square roots are precisely those that one 
would take on the Hilbert space TCl^ if the gravitational field was a background field in order to arrive 
at the annihilation and creation operator decomposition. 



(4.8) 



Let us now proceed to the details: 



Maxwell Hamiltonian 



The electromagnetic GauB constraint operator applied to cylindrical functions over 7 reads 

G^Q{A),= E ^(^)[ E E ^^]- (4.9) 

veVi'y) eeE(7);6(e)=i) eeE{'Y);f(e)=v 

Since we are working with gauge invariant functions, the Gauss constraint is identically satisfied. 
Using the notation Y\v) := Yej^v) we obtain the operator identity 

+ YjA^)] = J2lY'iv) - Y\v - I)] ^ {djY^){v) = 0, (4.10) 

where naturally the backward edge derivative has popped out (we used l^g-i = —Yf.) and Einstein's 
summation convention is implicit. 

Next consider the loops /J^^. ^2,0-3- It is easy to check that with euK = 1 for fixed I we have 

ln(^/3|i;.2,.3(^)) = ^lil'^(^ejK)) + HHeK(v'+J)) - HHej{v'+K)) ' HHeK{v')]^'=y+Z2^J+Z2^K- 

(4.11) 

Using the notation Ai{v) := [P"*" ■ ln(iJe7(i)))/'^]niod(27r) with P-^ defined in ( f4.16| ) we can rewrite 
(|4.11|) as (subject to the remark after ( [4.2|) ) 

HH,i^.^.,Jv)) = zaie^*^^(9+i^)^,^^^^,^^^ (4.12) 
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where naturally the forward edge derivative has appeared. Equation (f4.12|) is obviously gauge in- 
variant under Aj Aj + dfF. 

It is important to notice that forward and backward derivatives commute with each other, 
[df,dj] = for any I,J,a,a'. We can now introduce the one-particle Hilbert space Ti-l^^^ with 
inner product 

< F,F' >T^.^ = Y,FW)Fi{v) (4.13) 

V,I 

and one easily checks that {dfY = —dj"^ is the adjoint of the lattice derivative on Hl^ j^. 
It is convenient to introduce the lattice Laplacian 

{Af){v) := J2(97dtf){v) = + bj) + f{v - hj) - 2f{v)] (4.14) 

which is easily seen to be negative definite on i 

<F,AF>^^,=- Y.\^d-j Aj){v)\' (4.15) 
t^ey(7) i,J 

and invertible since does not contain zero modes by definition (they are not normalizable). 

With its help we may define the lattice transversal projector 

(P± ■ FUv) = Fj{v) - [di^djF']{v) =: P(t7),K,j)^'(^')- (4-16) 

v' ,v 

We then obtain the operator identities 

Y\v) = (Px ■ YYiv) and e''"" {d+ Ak){v) = e''"" {dj[P^A]K){v). (4.17) 

It is important to realize that the lattice metric Sjj is not a background structure, but is actually 
diffeomorphism invariant, it is the same for all cubic lattices and only depends on the topology of the 
lattice (which in our case is cubic). Therefore the index position of the index / is actually irrelevant, 
in particular, P± = P^. Cubic graphs are distinguished by the fact the same projector P± in 
maps to the space of solutions to the Gauss constraint djF^ = and to the gauge invariant piece of 
Fi under Fj ^ Fi + d+f. Notice that indeed P^ = P^ = p[, dijPtA 

v,v') = 25yy is a symmetric 

projector on i on the two physical degrees of freedom per lattice point as desired. It is remarkable 
that all the structure that comes with df can be constructed without any reference to the gravitational 
degrees of freedom). Of course, this is due to the fact that the exterior derivative of a one form and 
the divergence of a vector density are metric independent. 

Let us write the Hamiltonian Hm,'^ in our new notation. In order to simplify our life for the 
exploratory purposes of this paper we will replace the sum over the four loops corresponding to the 
choices cxi, a2 divided by four in ( [1.2|) by one loop corresponding to ui = (72 = 1- Likewise we replace 
the sum over the choices a divided by two by the term corresponding to a = 1. This just corresponds 
to the exploitation of the quantization ambiguity from which all the Hamiltonians constructed so far 
suffer anyway. Then ([4.2|) can be written in the compact form 



X 

~i ■• '2' 
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X {[-P^ . Y]\v)[P^ . Y]\v) + [e^^^9+(Px ■ A)r]{v)[e'''''dUP^ ■ A)^]{v)} (4.18) 
where Q{{v^r) = Q\^j^{v,r). Let 

Qij{v,r) := [Q{{v,r)]^Q'j{v,r) = -Q{{v , r)Q^j{v , r) . (4.19) 



Notice that while the operators Q\{v,r)^Q-^j{v^r) do not commute, in ([4.18|) only the the symmetric 

^t'kw,j) as operators on (^H\,^^ 



piece of Qij{v, |) survives in ( [4.18| ). We now define G^'i , , j-., G'^'j. . , j. as operators on 'H^ ^H]^. 



v'). 



by 

v',K 

[G^^-'^ij ® F]{A, {v, I)):=-J2 PiiiU<j/'''d;,MLM{v', ^)^)(A)6^'^^^(9+(P^ ■ P)p)( 

v',K 

We can then write the Maxwell Hamiltonian in the even more compact form 

< = ^[< y\ G^'-'y >ni,,, + < A\ G^^-n >^.^j (4.20) 

where the adjoint in ( [4.20| ) is with respect to Ti^^ . The following results are crucial. 



Theorem 4.1 

The operators G^^'^, G*^'^ are positive semidefinite and definite on the subspace ® '^1/^7 where 



a) 

The operator is positive definite on Ti^ ® TH.^ and thus H^^ is positive definite on Ti. 
Proof: 

i) 

Let ^ = E;. z^i)'^ ® F^, en^® nlj^^ be given. Then 



E E 1 1 E ■ ^^)'(^) [^K^, ] I \'n^ (4.21) 



where we used that P = P± commutes with edge derivatives. By the same manipulations we arrive 
at 
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= EEiiE^ME(^-^^><^^)'(^)[^>'^)^f]iiHf (4-22) 

V j fl I 

where {d"^ x F)^ = e^'^^djFx- The definiteness statement is clear by inspection. 



Let now ^^^^ = Zf^tpj^ (g) ^^^^ eTi^ (S) Hm^ be given. Then 



(1,1/ 

fl,U v,I,J 

+ < (P • a+ X (P ■ 5+ X J 



2 



+ EE II E^m[^?K^, ® [(^ ■ X (4.23) 

The definiteness statement follows from the fact that we are working on the space of gauge invariant 

functions (solutions to the Gauss constraint). 

□ 

It follows from this theorem that the operators G*^'^, G*^'^ have an inverse on the subspace Ti^^Ti]^^. 
This fact is vital in order to arrive at the creation and annihilation operator decomposition of the 
Einstein-Maxwell Hamiltonian. We state here without proof that the positivity property of H!^^ 

holds on every graph so that it extends to the whole Hamiltonian H^^ . The latter statement follows 
from the fact that by construction the consistently defined Hamiltonian (|3.40|) preserves the space of 
gravitational spin-network functions over a graph due to the projection operators on both sides in 
(|3.39|) . (Actually this is already true for the operator defined on spin-network functions because our 
operators do not change the graph on which a state depends). Therefore, if \l/ = ip^^ ^ip^l + ip^^ 
we have 



< VI/, H^'^ >n,.„=< < ® ® < > + < < ® i^Z^KK ® < > (4-24) 

due to orthonormality of the spin network functions and both terms are separately positive. 



Klein-Gordon Hamiltonian 

Here the one-particle Hilbert space 'Hj^Q.^ is equipped with the inner product 

(^'^')«u.= E n^W)- (4-25) 

v€V{-y) 

Let us also simplify the Klein- Gordon Hamiltonian ( |4.3|) by replacing the average over the eight right 
oriented triples of edges by a single one corresponding to ai = a2 = as = 1 

Hkc-t = — ^ii^'^P 2^ ^ 
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DeV(7) 



(4.26) 



One can read off from ( [4. 261) the operators G^^,'^, G^^'^ on giving rise to operators G^^'^, G^^''^ 
on ® '^KG,'y which allows us to write ([4.26|) in the compact form 



KG 



(4.27) 



where the the dagger is with respect to Ti:^ and we used the notation (j){x) = \n(U{x))/i. Of course 
there is no projection operator involved in this case since in this paper we deal with neutral scalar 
fields only. 



A theorem analogous to theorem |4.1| can be proved in this case as well and will be left to the 
ambitious reader. 



Fermion Hamiltonian 



Finally, let us also simplify ( |4.4] ) by replacing the average over the eight triples of edges to the 
one corresponding to ai = (T2 = as = 1, that is 



H 



D,-y 



mp 
'2K 



J2 [OB{v')ei{v)~W)e'l{v)] 



v,v'£V(-y) 



X {{^ijk^^'^^Q\{v, ^)Qj(w, ]^)[r^{hK{v)5^>^^+K - S^',v)\ab} 



(4.28) 



where hi{v) = hej{v)- One can read off from ( |4.28| ) the operator G^^^^ j-^, ^^-j on giving rise to 

the operator G^ on ® '^z?,^- Here /x = 1, 2 where 9\{v) = Oa{v) and 0\{v) = 0'^{v) and Ti-j^.^ is 
equipped with the inner product 



(4.29) 



This allows us to write ( ^.281 ) in the compact form 



-D,7 



(4.30) 
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where the the dagger is with respect to Ti^l . Of course there is no projection operator involved in 
this case since in this paper we deal with neutral spinor fields only. 

A theorem analogous to theorem [4.1| has not been proved in this case and is fortunately not 
necessary in order to arrive at creation and annihilation operators. See the remark at the end of 
subsection O . 



4.3 Fundamental Fock States and Normal Ordering 

We notice that both bosonic Hamiltonians have the structure 

^7 = ^ E + (4.31) 

where pi{v), qi{v) are operator valued distributions on 'Hmauer,'-/ with canonical commutation relations 
[pi{v),qi{v')] = ih6vy6iii and P{{v,l),{v',l')),Q{{v,l),{v',l')) define positive definite operators on 
T^geo,^ ® H^- In particular, they are symmetric there, that is, P{{v,l), {v' ,l')y = P{{v' ,1'), {v,l)) 
where the dagger is with respect to TCgeo,'y It follows that for a state ip G 7igeo,7 expectation 
value 

P^iv, I), iv', I')) :=< i;, Piiv, I), K, 0)^ >H,e... (4.32) 

defines a positive definite, Hermitian matrix on 7i^. Of course, we are being here rather Cavalier 
concerning domain questions and self-adjoint extensions (one possible choice is the Friedrichs exten- 
sion) but a more detailed analysis would go beyond the exploratory purposes of this paper. (If S 
is compact then 7 is a finite graph and all operators are bounded, although then we need to talk 
about boundary conditions). With these cautionary remarks out of the way we then can state the 
following theorem. 

Theorem 4.2. 

Suppose that the operators P{{v,l),{v'J')),Q{{v,l),{v',l')) form an Abelian subalgebra of C(Hgeo,'y) 
and that they are self- adjoint for each pair {v,l), {v',l'). Then there exists a unitary operator U on 
'Hgeo,-y ® Ti.^ such that 

^ [< p, Pp >ni^ + <q,Qq >ni^] =< z\ ujz >ni^ 
h{v) = -^Y.ia{ivJ),{v'J'))q,iv')-fbiiv,l),{v',n)piiv')] 

a = lib 



b = {d-Y 

oj = (a 



-^^^Qd^^ (4.33) 



where the square roots and inverses are with respect to 'Hgeo,-) ®'H^ while transposition is with respect 
to Til^. Moreover, if we set 

i=[5^[a((t;,/),(t;',r))'?r(^^')-^&((^,0,(^',O)P/'(^')], (4-34) 
v2 ,„ 



v'l 
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then Q(v),c|(f) satisfy the canonical commutation relations where the adjoint is with respect to 

a) 

Suppose that we are given real valued and symmetric operators P,Q on Til^. Then there exists a real 
valued, unitary operator U on such that 

^ [< Pp >^i + < g, Qg >^i ] =< z, ujz >-j^i 

z = -^[aq-ihp] 
a = UD 



h = {a-Y 

oj = {a~YQa-^ (4.35) 
where the square roots, transposes and inverses are with respect to to 7i\. Moreover, if we set 



ci{v) 



i= 5^[a((t;, /), {v\ l'))qv{v') - ih{{v, I), {v', I'Miv')] (4.36) 

v'V 



then ci{v),cl{v) satisfy the canonical commutation relations where the adjoint is with respect to 

T~^matter,'y ■ 

The proof is straightforward and is omitted. Notice that unitarity means that (f/t)^ = [/-I where 
the dagger is with respect to 'Hgeo^-y while U'^ = U~^. The second assumption in i) is not immediately 
satisfied for the operators P, Q in ( [4.18|) and (|4.26| ) as they stand because in the form displayed they 
are only symmetric, (-P^)"^ = P, and similar for Q. However, since, with respect to the Ti}^ degrees 
of freedom they act on matter operators of the form p ® p and since the p commute among each 
other, we may without loss of generality assume that P^ = P so that P = P'^ . The first assump- 
tion in i), however, is violated for the geometrical operators P, Q of (|4.18|) and ( [4. 361) , they do not 
commute on general states. On generic states, however, they do commute. This non- commutative 

geometry will lead to further quantum corrections for what follows, meaning that the operators c, c 
satisfy canonical commutation relations on generic states only but not exactly. We will not discuss 
these effects in this paper and from now on assume that P{{v, I), (f I')), Q{{v, I), {v', I')) generate an 
Abelian operator algebra on Hgeo,'y 



With these cautionary remarks out of the way, theorem [4.2| suggests to choose a different order- 
ing for the operator (|4.31|) , namely the normal ordered form 

^7 = 5Z ^ii^y'^{v,i),(v',i')Ci'iv'). (4.37) 

iv,l),iv',l') 

When comparing ( |4.37|) with ( [4.31| ) one finds out that they differ by a purely gravitational operator 
which is the quantization of the usual, IR divergent, normal ordering constant in flat space. This 
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can be avoided as follows: The clean way to arrive at the form ( [4.37 ) from first principles is to write 
the classical expression 

^ = ljj''^j^ A[^((^, 0, (y, l')Pi{x)pv{y) + g((x, /), (2/, l')qi{x)qv{y)] (4.38) 
whose quantization gives rise to ( [4.31|) , first classically in the form 



H= d-'x d-'yci{x)uj{{x,l),{y,l'))c,{y) (4.39) 

and then to quantize it. However, in order to do that we would need, for instance, the explicit 
expression of the functions uj{{x, I), {y, I') in terms of the elementary gravitational degrees of freedom 
Ai^, Ej which is unknown. Thus, our procedure to first quantize (|4.31| ) and then to normal order it 
should be considered as the "poor man's way" of quantizing ([4.39|) directly which would not lead to 
a normal ordering operator. 

We have judiciously chosen the double hat notation for the operator c in order to indicate that 
it involves the anticipated mixture of gravitational and matter quantum operators. 

We suggest that the operators ci{v) play the role of the fully geometry — matter coupled 
system that is normally played by the matter annihilation operators on a given back- 
ground geometry. 

In order to justify this, we should now construct coherent states of 'Hgeo,'y ® 'Hmattera that are eigen- 
states of c and from them a vacuum state and n-particle states. (Notice that such states are auto- 
matically embedded in the full kinematical Hilbert space). We will choose the complexifier method 
p8| in order to do that. 

The idea of a complexifier is to find an operator (7, which in our case will depend on both 
gravitational and matter degrees of freedom, such that 



e 



(As before, we are working on Ti^ = 'Hgeo,-y ®'Hmatter,'y for each 7 separately). Comparing with (|4.33| ) 
we find the unique solution 



C = Cg,, + \ J2 Pi{v)D-Mv,l),{v',n)p,{v') (4.41) 



2 

{v,l),iv',l') 

where Cgeo is a positive definite operator constructed from the gravitational electrical degrees of 
freedom only according to the guidelines of so that [C — Cgeo, Cgeo] = 0. 

The complexifier coherent state machinery can now be applied and we arrive at the following 
coherent states: Let m be points in the full phase space of gravitational and matter degrees of 



freedom. Let C be the classical limit of ( 4.14| ) which we know in principle exactly in terms of 



Q, P which were classically given in terms of the gravitational three metric and partial derivative 
operators. Compute ^^^(m) = [e~'^'~'^c A^]{m) and2;m(^) = [e~*^'^cg](m) where is the gravitational 
connection, q = or q = (p are the Maxwell connection and Klein Gordon scalar field respectively, 
Xc is the Hamiltonian vector field of C and C denotes the Lie derivative. Both functions are functions 
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of both matter and geometry degrees of freedom. Let 5/i';7®'^//';7 be the 6 distribution with respect to 
the uniform measures of the L2 spaces that define 7ifcm,7 denote by h', H' the set of gravitational 
and matter (point) holonomies along the edges of 7. Then, for instance for Maxwell matter 

V'm;7 '■= {e'^^^' [5h';7 ® ^ H' h' -^h{zE(m),H' ^H{zM(m)) (4.42) 



define coherent states on Hkin,^ as shown in ||38[. Here h{zE{rn)) denotes the set of gravitational 
holonomies h{A^) along the edges of 7 where the real connection A is replaced by the complex 
connection ZEim) (analytical extension) and similar for H{zM{fn)). One of the nice features of the 
coherent states (|4.42| ) is the fact that they are eigenstates of the operators ci{v) with eigenvalue 
Q(f)[m] as one can explicitly check (using the fact that d commutes with D). 

As an example, consider the case of photons propagating on fluctuations around fiat (i.e. empty) 
space. Then we have := {Ai^,Ej) = (0,5°) and := {Aa,E"-) = (0,0) so that for an edge 
e G -^(7) we have He{zM{jnP)) = 1- We have qiiy) = qi{v) = ln(iJe/(»;))/^ so that Q(f)[m*^] = 0. Thus, 
the operators ci{v) annihilate the vacuum state fi^ := ^mO-rf 7- Moreover, since [^(f ), Q/(f ')] = 
f^v,v'5i,v we are able, in principle, to construct a symmetric Fock space J-'^{Hl^) where the n— particle 
states are defined by 



1/1, . . . , /n) := aVi) • • • c{F^)n, where 1\f) := J] Fi{v)^,{v). (4.43) 

v,l 



A precise map between (|4.43|) and the usual photon states on fiat Minkowski space can be given for 



the fundamental states ( |4.43| ) as well but we postpone this to the next section. 



Notice that due to commutativity of different matter types we can add the operators C in (|4.41| ) 



for different matter types in order to arrive at simultaneous coherent and Fock states for all matter 
types! 

4.4 Approximate Fock States and Semiclassical States 



It is clear that the program sketched in section 4.3 cannot be carried out with present mathematical 



technology because we are not really able to construct the operators a, b, to which require precise 
knowledge of the spectrum of these operators on 'Hgeo,'y ® 'H}^. Thus, in order to to proceed we have 
to do something much more moderate. As a first approximation we consider states which do not 
mix matter and geometry degrees of freedom in the way ( [4.42] ) did but rather will look for Fock 



states of the form ipmEi-y ® V'matter;7 where ipms-n ^ gravitational coherent state peaked at the point 
G A^£; in the purely gravitational phase space, for instance the ones constructed in [1^, |ll|]. 
These states are generated by the piece Cgeo of the complexifier of (|4.41| ) by applying its exponential 
to the 6 distribution on the gravitational Hilbert space alone. Let {T„} be a complete ortho normal 
basis of states in 7-^mattcr,7- Using the overcompleteness of the just mentioned coherent states, we 
can write the matrix elements of H.y, given in the normal ordered form of (|4.37 ) as (interchange of 
summation and integration must be justified by closer analysis) 

< i^ms-n ® ''/'matter; H^t/j^^-j ''/'matter > 

ruE-n ® '/'matter; ''PmEl'y ® ^ ^ 

{v,l),{v',l') ''^E ''^E n,n' 
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® Matter > • (4-44) 



Here z/^ is the Hall measure |^ generalized to graphs in [q5| and Ai'^^ is Ai restricted to the graph 
7 as defined in p6| . 

Let us introduce the real valued operators on Til^ defined by 

a"'^{{v,l),{v\l') ■=< ^^,,a{{v,l),{v',l'))^mE >h^ (4-45) 

and similar for a), h. Consider also the operators on 7imatter,7 defined by 



v'V 



The coherent states 'ipmE-^ ^'^^ sharply peaked in M."' which implies that up to h corrections 

® ^matter >= ^u, (^ij, m'^) < T/^^atter, c'"^Catter > (4-47) 

and similar for a). We conclude that up to h corrections 

< '^ruE-a ® "^matter, H^tpmE;j ® V^matter > 
= E <C"(^)^matten<!),(.,.)ar(^')^Ltter>- (4-48) 

Now, again using that the ipmE-n have very strong semiclassical properties, it is possible to show that, 
up to h corrections, the operators a™^ , 6"*^ , u;'"^ can be computed by first calculating the expectation 
values of the operators P, Q on ®H}y, which we know explicitly, to arrive at operators P'^e ^ q^e 
on Tt}^ and then to plug those into the formulas (|4.35| ). Thus, we have arrived at a "poor man's 
version" of an annihilation and creation operator decomposition of which approximates the exact 
version but which still takes the fiuctuating nature of the quantum geometry into account through 
the uncertainties encoded into the states 'ipmEi-y- 

It is now clear how we arrive at approximate Fock states. Instead of (^4.41|) we consider the 
operator on 7imatter;7 defined by 

C"'^ = l E Piiv){D"^n-Mv,l),iv',l'))p,iv'). (4.49) 

(v,l),{v',l') 

This complexifier now generates coherent states on 7^matter,7 in analogy to ( [4. 421 ), e.g. for Maxwell 

matter, by 

■= {e-^^'"'5H'-a)H'-.Hi.-HmM)). (4-50) 

where mj\/ is a point in the matter phase space, z^'^irriM) = {e~^^^c"^E q^(^rnM) and C'^'^{mM) is 
the classical limit of C*™^ on the matter phase space. Choosing the points m%j appropriate for 
vacuum will now produce a vacuum state := i'^o for the operators ff^E which by construction 
(theorem [4.2|) satisfy canonical commutation relations exactly. Summarizing, with 

^ruE-,^ = (e-^-°/%0/.'-M4{-^))' ^'e(^i^) = e'''''^^-°A^ (4.51) 
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we arrive at approximate ri— particle states 

...,E,y = ^m,,-y ® &\F,) . . . c\Fr,)Q'^, (4.52) 

where c^{F) = J2v i Fii'^)^K'^) associated Fock space. 

Let us conclude this section with some remarks: 

For the lepton sector things are much easier because the operators 0'^{x), (^^4(3;))^ already satisfy 
canonical commutation relations and all momenta are ordered to the left in H^^. A suitable vacuum 
state annihilated by all the ^^(f) is given (up to normalization) by 

^^(0) = n ^a(^)- (4.53) 

v,A,fj, 

2) 

In order to relate, say the states ( [4.52| ) to the usual Fock states on Minkowski space with Minkowski 



vacuum Qp and smeared creation operators Ci?(/) = / (fxfL{x){cpy{x) with different labels L we 
need the particulars of the expectation values of the gravitational operators. Basically, the discrete 
sum involved in the definition of c{F) is a Riemann sum approximation to the integral involved in 
the definition of cp{f) which becomes exact in the limit that the lattice spacing e vanishes where 
F/{v) := e"X/"(f )/i(i;) for some power of e and some matrices X[" which depend on the choice of 
the gravitational coherent states. One can then establish a map between the usual Fock states and 
our graph dependent ones by 

4(/i) . . . cUfn)np ^ c\Ff^) . . . c\Ff-)n, (4.54) 

which becomes an isometry in the limit e ^ 0! In more detail, let us consider the simple example of 
a regular cubic lattice in R^. Discarding fluctuation effects from the gravitational field we arrive at 
the dimensionfree graph annihilation operators for Maxwell theory given by 

c'^{v) = 'y=A;i, - 1 'y=A;''E^](t;), (4.55) 

y la 

where = . On the other hand, for the usual Fock representation we get the annihilators 

of dimension cm~^/^ given by 



"la 



^[ ^W^Aa - I 't/^ '^T(x). (4.56) 



Using dimensionfree transversal fields F\v), djF^ = on the polymer side and transversal fields 
/"(f), daf^ = of dimension cm~^/^ we arrive at smeared, dimensionfree creation operators of the 
form 

c\{F) = J2 F\v)c'^{v) and c^/) = / d\ r{x)cUx). (4.57) 

Using E\v) ^ e^^a-E^, Aj{v) ^ e5jAa{x), A^{v) e'^A{v) we see that c!^{v) ^ e^^'^5lcp{v) so that 
our desired map is given by 

(F^Yiv) = 5iiv)e'/'nv), (4.58) 
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where the two factors of e^/^ combine to the Riemann sum approximation of the Lebesgue measure 

Finally we mention that the way, approximate n-particle states were obtained above, bears some 
similarity to the treatments of gravitons in |T2|, [1^. In both semiclassical state for the 

gravitational field is used to obtain a classical background geometry. In |jT2|, Q a weave state is used 
for this purpose, here we have employed the coherent state ipmE- other respects, the treatments 
differ, however. To define the notion of gravitons, a split of the gravitational field in a dynamical and 
a background part is necessary whereas nothing of this sort is required for the treatment of matter 
fields. 



5 Towards Dispersion Relations 

Dispersion relations are the relations between the frequency uj and the wave vector k of waves of a 
field of some sort, traveling in vacuum or through some medium. In quantum mechanical systems, 
the dispersion relation is the the relation between the momentum and the energy of particles. The 
form of the dispersion relations appearing in fundamental physics is dictated by Lorenz invariance. 
Since this invariance is likely to be broken in quantum gravity, modification of dispersion relations 
is conjectured to be an observable effect of quantum gravity. In this section we would like to explain 
why QGR indeed leads to modified dispersion relations, and how one might proceed in a calculation 
of these modifications. 

There are at least two mechanisms by which modified dispersion relations arise in the context of 
QGR, and it is important to keep them apart. Let us start to discuss the first one by considering an 
analogous effect in another branch of physics: 

A prime example coming to mind when thinking about modified dispersion relations is the prop- 
agation of light in materials. The mechanism which causes these modification is roughly as follows: 
The electromagnetic field of the in-falling wave acts on the charges in the material, they are acceler- 
ated and in turn create electromagnetic fields. These fields interfere with the in-falling ones, the net 
effect of this is a wave with modified phase and therefore, a phase velocity differing from the one in 
vacuum. The precise relation between the force acting on the charges and the fields induced by them 
depends on the properties of the material and also on the frequency of the wave, and thus gives rise 
to a frequency dependent phase velocity and, hence, a nontrivial dispersion relation. Under some 
simplifying assumptions, this relation looks as follows: 

where ujq, k, and p are properties of the material. As is to be expected, if the energy of the in-falling 
wave is very low compared with the binding energies (~ ujo) of the charges, the frequency dependence 
of the phase velocity will also be very small. 

In QGR, modified dispersion relations can be expected from the interplay between matter and 
quantum gravity by an analogous mechanism: The propagating matter wave causes changes in the 
local geometry, which in turn affect the propagation of the wave. Again, if the energy of the wave is 
very small, so will be the modification of the dispersion relation as compared to the standard one. 

In order to honestly account for this back-reaction mechanism one would have to do a first 
principle calculation that involves solving the combined matter - geometry Hamiltonian constraint. 
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Technically, we are not yet in the position to do that. However, as a first approximation we can take 
care of the reaction of the geometry to the matter fields by using coherent states '^ms^ ^matter in the 
constructions of section ^ which are peaked at a classical configuration which is a solution of the field 
equations of the combined gravity-matter system. This should be seen in analogy to our remarks in 
section 2 where QED corrections are computed with coherent states for free Maxwell theory instead 
of the full QED Hamiltonian which neglects the back-reaction from the fermions. 

There is, however, a second source of modifications to the dispersion relations: The inherent 
discreteness of geometry found in QGR. This effect has nothing to do with back-reaction of the 
geometry on the matter and it is the contribution of this effect to the dispersion relations that can 
be studied with more confidence. This is what we will discuss in the rest of this section and in our 
companion paper 

Let us again start by briefly reviewing an analogous phenomenon from a different branch of 
physics, the propagation of lattice vibrations ("sound") in crystals. As an example, consider an 
extremely simple model, a one dimensional chain of atoms. We assume that all atoms have the same 
mass m and that each of them acts on its two neighbors with an attractive force proportional to the 
mutual distance. If we denote by e the interatomic distance in the equilibrium situation, by q{z) the 
displacement of atom z from its equilibrium position ez and set p{z) = q{z), the Hamiltonian for the 
system reads 

H = IJ2-p'(')+^ (g(^ + 1) - q{z)f . (5.1) 

The corresponding equations of motion are simple, a complete set of solutions is given by 

2K 

q(t, z) = expi (ezk — Lj(k)t) , with u'^ik) = (1 — cos fee). (5.2) 

m 

As the solutions are straightforward analogs of plane waves in the continuum, uj'^{k) is readily inter- 
preted as the dispersion relation for the system. We see that it contains the "linear" term propor- 
tional to fc^ expected for sound waves in the continuum, as well as higher order corrections due to 
the discreteness of the lattice. 

Let us reconsider (|5.1| ): The fact that the q{z) are displacements of atoms is not explicitly visible. 
H could as well be the Hamiltonian of a field q with a certain form of potential, propagating on a 
regular lattice. Having made that observation, we are already very close to the model just described. 
Upon choosing a semiclassical state \l/ for the gravitational field, the bosonic Hamiltonians of section 
^ are of the form 

H^ = \ E p/(^)i^*((^,o,(^',n)p.(t;')+g/(^)Q*((^,o,(^',n)^/'K), (5.3) 

where the expectation values 

p^{{v, 0, {v\ I')) = {Piiv, 0, iv\ n)h Q^iiv, 0, iv', n) = {Qiiv. 0, n))^- 

contain an imprint of the fluctuations of the gravitational field. can in principle be taken to be 
a coherent state for the gravitational field peaked at an arbitrary point of the classical phase space. 
However, since we are interested in dispersion relations, a notion that by definition describes the 
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propagation of fields in fiat space, we will restrict considerations to the case of GCS approximating 
flat Euclidean space (denoted by ^fflat in the following) 

There are however two fundamental differences between (|5.1| ) and (|5.3|) , and we will discuss them 
in succession. The first difference is that ( ^.31 ) is a Hamiltonian for a quantum field whereas the 
former is purely classical. Note however, that the Hamiltonians of section ^ are normal ordered. 
Thus, the expectation value of these Hamiltonians in a coherent state peaked at a specific classical 
field configuration (defined in sections [4.3| , [4.4|) will yield precisely its classical value. Moreover, 



expectation values for the matter quantum fields will exactly equal the the corresponding classical 
values at all times. Therefore, in discussing the dispersion relations, we will assume the matter 
quantum fields to be in a coherent state and can effectively work with the classical fields p, q. This 
will be a very good approximation in processes such as light propagation due to decoherence of 
matter quantum effects in large ensembles of photons and because higher loop corrections of QED 
are not Poincare violating (which is what we are interested in here). 

The second and more important difference between ( |5.1|) and ( |5.3| ) lies in the following: In (|5.1|) , 
the coefficients of the fields do not depend on the vertex. This is the reason why one can explicitly 
calculate solutions to the equations of motion. In contrast to that, P{{v,l), (■ ,/')) will in general 
depend on v, even if the state \l/fiat employed to compute the gravity expectation values is a good 
semiclassical state. As a result, the field equations will be complicated and, most important for us, 
not have "plane wave" solutions 

q^(t, v) = expi{kx{v) — ut) (5.4) 



any more. Hence if we would Fourier decompose solutions of the field equations with respect to ( |5.4|) , 
the support of the resulting functions will not be confined by a dispersion relation to some line in 
the u)-\k\ plane, anymore. 

However, for a good semiclassical state, symmetry, which is absent due to the vertex dependence 
of the coefficients, will be approximately restored on a large length scale. For example, if the 
vertex dependent coefficients would be averaged over large enough regions of S the average would 
be independent of the specific choice of the region. Therefore, for long wavelength, plane waves ( |5.4| ) 
should at least be approximate solutions to the field equations. The following scenario is conceivable: 
Although there is no exact dispersion relation, the support of the Fourier transform of a solution 
might be confined to some region in the uo-\k\ plane, or the Fourier transform has at least to be 
peaked there. This region should get more and more narrow for longer wavelength, leading to an 
ordinary dispersion relation in the limit (see figure |^). We have to note, however that even if this 
is true, there is no guarantee that a dispersion relation with corrections to the linear term makes 
sense as an approximate description for long wavelength. We tried to visualize this in figure ^ 
So, to conclude, it is very plausible that a nonlinear dispersion relation will turn out to be a good 
approximate description of the physical contents of ( |5.3| ) for long wavelength in this sense. But issues 
such as the one depicted in figure ^ definitely merit further studies. 

Let us now turn to the practical question of how a nonlinear dispersion relation can actually be 
computed from (|5.3|). A simple model that displays some of the complications due to the vertex 



dependence of the coefficients in (|5.3|), but can nevertheless be treated with analytical methods, can 



•^Also, when considering application to situations such as the 7-ray burst effect, the curvature radius is always huge 
compared to Planck length and does therefore not lead to any new quantum effects but just to classical redshifts which 
can easily be accounted for. 
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Figure 2: Fourier transform of EOM. The support of solutions has to he in the shaded region, 
a) b) 





Figure 3: Can higher order corrections to the dispersion relation be given? a) Yes, approximately, 
b) No, there is no meaningful notion of dispersion relation beyond linear order 



be obtained from (|5.1| ) in the following way: Upon setting m = 1 and K = the Hamiltonian 
(|5.1|) can be interpreted as that of a scalar field propagating on a one dimensional lattice with lattice 
spacing I. Now we partly remove the assumption of a constant lattice spacing by replacing I hj Iz - 
the distance between the lattice point labelled by z and the one labelled hj z + 1 - but still assume 
periodicity on a large scale. 



h = l 



N- 



for all z G Z 



for some iV G N. It turns out that the dispersion relation for this system has several branches and 
the small k behavior of the acoustic branch is given by 



12 ((F)) 



\kf + 0{\kf) 



(5.5) 



Here, 



Af-l 
n=0 



N-1 



{j -i)[N- (j - i)l and ((/)) = ]^Y.^n etc. 



n=0 



As we we have indicated by means of notation, it is instructive to view the /j, i = 0, . . . , — 1 
as independent random variables. The moments of the corresponding distribution determine the 
dispersion relation ( |5.5|) . In each order in fc, corrections are present as compared to the case of 
constant lattice spacing. Note in particular that the phase velocity \im.k_,Q uj {k) /k can be smaller 
than one. 
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It is remarkable how subtle the dependence of the dispersion relation on the distribution of the 
lengths is already in this simple model. Although we do not have a proof, it is plausible that quali- 
tatively, the above formula extends to higher dimensions, where analytical proofs get much harder. 
For more information about the model as well as a proof of ( p.5| ) we refer the reader to or 



for a beautiful numerical analysis of similar models in two dimensions to p9 



The discussion of the one dimensional model given above shows how complicated an exact analysis 
of the equations of motion is already in simple cases. Since the models 

^*fl. = ^ J] pK^)i^*fl.((^,0,K,O)pK^') + ?/(^)Q>.fl.(^,0,(^',O)*'K), (5.6) 

v,v',l,l' 

that one obtains from QGR are more complicated, it is useful to explore a less precise but easier 
route towards dispersion relations. The idea which we would like to advocate is to replace (|5.6|) by 
a simpler Hamiltonian which 



is a good approximation of ( |5.6| ) for slowly varying q and p and 
is simple enough such that the EOM can be solved exactly. 



This idea underlies also the works and and, at a rather simple level, is the basis for the 
recovery of continuum elasticity theory from the atomic description in solid state physics (see for 
example pO[ ). 

We will now propose a replacement for ( ^.6|) fulfilling the above requirements. In the long wave- 
length regime, we can revert to the continuum picture, i.e. replace the lattice fields qi{v),pi{v) by 
qi{x{v)) , pi{x{v)) , where now qi,pi : S ^ C. 

We should however take care that information about the lattice is at least partially encoded in 
the continuum theory. To this end, we Taylor expand the (now continuum) fields in the Hamiltonian 
up to a certain order. For example, we would make the replacement 



qi{x{v))Q^i{v,l),{v',l'))q,{x{v')) 

1 
2 



qi{xiv))QAivJ)Av'J')) b\v,v')d,qiixiv)) + -b'iv,v')biiv,v')did,qiixiv)) 



where U{v,v') = x{v') — x{v). Even if we terminate the Taylor expansion after a few terms, the 
resulting Hamiltonian will be an excellent approximation to the original one, provided p and q 
change only very little from vertex to vertex, our standing assumption in the whole procedure. 

Now we will eliminate the spacial dependence of the coefficients, which was the main difficulty 
in dealing with the original Hamiltonian (|5.6| ), by replacing them by their averages over all vertices 
of the graph. This can be justified as follows: If the fields p, q are varying considerably only on 
length scales much larger then some macroscopic scale L, it is a very good approximation to replace 
the vertex dependent coefficients in the Hamiltonian by their averages over the vertices in regions 
of dimension L^. On the other hand, as we have said before, a good semiclassical \l/fiat state will 
insure that the system described by ( p. 6] ) has the symmetries of fiat space at least at large distances. 
One way to state this more precisely is that the average of the coefficients appearing in ( |5.6|) over 
vertices in regions with characteristic dimension or larger, is independent of the region to a good 
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approximation. Therefore we can indeed replace the vertex dependent coefficients by their averages 
over all vertices. 

Let us again give an example for a typical term in the Hamiltonian: 

J2QiivJ),{v',l'W'{v,v') ...b^'^{v,v') =: F'''-^"{v,l,l') > {{F""{-,l,l'))) 

v' 

where we have introduced the graph average 

/,/'))) = (5.7) 

V 

N being the number of vertices of the graph \l/flat is based on. In case we are dealing with an infinite 
number of vertices, the definition (|5.71 ) has to be replaced by the limit of averages over finite but 
larger and larger numbers of points. 

Finally we can replace the sum over vertices of ( ^.4|) by an integral. We thus end up with a 
Hamiltonian for a continuum field theory on E and the coefficients of the fields being constant. 
Therefore the equations of motion of the theory admit plane waves as solutions, and their dispersion 
relation can be computed and discussed. This dispersion relation should describe the physical content 
of ( |5.6| ) for low energies (large wavelength). 

To justify our procedure, let us point out again, that it uses both assumptions (large wavelength 
- homogeneity and isotropy of the state on large scales) that seem essential from physical consider- 
ations, to recover a dispersion relation from ( ^.61) , enter in a transparent way. Also, if we apply the 



procedure outlined above to the simple regular lattice system ( |5.1| ), we recover, order by order, the 
nonlinear dispersion relation (|5.2| ). Thus, at least in this example, the simplified continuum theory 
still captures the information about the lattice to any desired order of accuracy. 

In the companion paper we will elaborate on the procedure described above and apply it to 
derive approximate dispersion relations from the Hamiltonians constructed in this paper, evaluated 
in the gauge theory coherent states of H. 



6 Summary and Outlook 

The goal of the present work was to begin investigations of the structure and semiclassical limit of 
the theory obtained by coupling matter fields to QGR. A basic assumption that we made was that 
the complicated dynamics of a full theory could be approximated by treating the matter parts in the 
Hamilton constraint of the full theory as Hamiltonians generating the matter dynamics and by the 
use of semiclassical states in the gravitational sector. 

Using this assumption we obtained the following results: 

1. We have proposed quantum theories of scalar, electromagnetic and fermionic fields coupled to 
QGR. The dynamics of these theories is generated by a Hamiltonian in the same way as in 
ordinary QFT. Consequently we were able to identify approximate n-particle states which cor- 
respond to the usual Fock states for matter fields propagating on classical geometries. In other 
respects, the theories are very different from ordinary QFT, thus reflecting basic properties of 
QGR: 
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• The basic excitations of the gravitational field in QGR are concentrated on graphs. The 
requirement of diffeomorphism invariance forces the matter degrees of freedom to be con- 
fined to the same graph as the gravitational field. The matter fields are therefore bound 
to become quantum fields propagating on a discrete structure. 

• In ordinary QFT, the background metric enters the definition of the ground state and 
the commutation relations of the fields. In QGR on the other hand, the geometry is a 
dynamical variable, represented by suitable operators. A QFT coupled to QGR therefore 
has to contain these operators in its very definition. This is refiected in the theories of 
section ^ by the fact that their annihilation and creation operators act on both, the one 
particle Hilbert space of the matter fields and the Hilbert space of the geometry. 

We also showed how a "QFT on curved space-time limit" can be obtained from this theory, 
using a semiclassical state of the gravitational field. 

2. We have discussed how modified dispersion relations for the matter fields arise in the context 
of QGR and motivated a method for computing them from the (partial) expectation values of 
the quantum matter Hamiltonians in a semiclassical state. 

Certainly, the present work can only be regarded as a first step towards a better understanding of 
the interaction of matter and quantum gravity. In future work, the assumptions that have been used 
should be removed, or their validity confirmed. 

On the other hand, application of the results of the present work can be envisioned. For example, 
it will be very interesting to see, whether the methods used in the present work can also be applied 
to investigate how gravitons arise in the semiclassical regime of QGR. This will be the topic of [1^ 



As another application, the companion paper contains a calculation of corrections to the standard 
dispersion relations for the scalar and the electromagnetic field due to QGR. 

To summarize, the interaction of quantum matter and quantum gravity is a fascinating but, alas, 
very comphcated topic, of which a good understanding still has to be gained. We hope that the 
present work illuminates the difficulties encountered in this endeavor and also contains some first, 
albeit small, steps towards its completion. 
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Appendix: Kinematical vs. Dynamical Coherent States: A 
Simple Example 

In systems with constraints linear in the basic variables, the expectation values of Dirac observables in 
a coherent state in the kinematical Hilbert space equal those in a dynamical coherent state, provided 
that both states are chosen to be peaked around the same point in the constrained phase space. This 
does not hold true anymore for systems with nonlinear constraints. One expects, however that the 
discrepancies between the expectation values on the kinematical and on the dynamical level will at 
least be small. In this appendix, we demonstrate that this is true for a simple quantum mechanical 
model system with a nonlinear constraint: A system of two coupled harmonic oscillators. 
Let the Hamiltonian of the harmonic oscillator be given as 

rj ^ f ^ 2 2 

H — - i h mu q 

2 \m 

It is well known that it can be quantized in terms of annihilation and creation operators a,a^, 
[a, a^] = 1 on the Fock space Ti over C. a is the quantization of the classical quantity 



[rmJ ( p \ 



A basis of H is given by the eigenvectors of the number operator N = a^a and will be denoted by 
\n). The coherent states for the harmonic oscillator are defined as 



/ 1^21 \ °o 



n=0 

A system of two harmonic oscillators can be quantized on the tensor product 7i^^^ = H ^Tt, the 
annihilation operators of the respective oscillators are given by 

Oi = a 1, 02 = 1 (8) a, 

and similarly for the number operators A^i, A^2- Analogously we have 

1^1,712) = |ni) (g) 1^2), \zi, Z2) = \zi) ^ \z2) , ^1,712 e No, -si,^2 e C 

and these vectors form dense subsets in H^™. 

Let us now impose the constraint C = Ni — N2 forcing the energies of the two oscillators to be 
equal. The kinematical phase space can be labeled by (-21,-22) G C^, the physical phase space by 
2; e C, where the embedding of the latter in the former is given by 

\z\ = |-2l| = I-22I , Z \z\ — -21-22- (A-l) 

The quantization of the constraint is simply 

a = TVi - N2, (A.2) 
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the physical subspace of H is given by 



n^^y^ = span{|n,n), n G Nq}. 
On this subspace, a new annihilation operator can be defined by 



on span{|n, n), n G N} and a^|0,0) = 0. It fulfills 
coherent states as 



CLka. Qj, 



1. Therfore we can define physical 



\z)(^ = exp 



n=0 



These are to be compared with the kinematical coherent states 1^1,^2); bearing in mind the identifi- 
cation ( |A.1|) . 



Let us consider the expectation values of the Dirac observables 0102 and A^^i and more complicated 
ones constructed from them. We start with A^^i: Clearly 



(A^i)|^i,^2) = \zi\ = \z\ . 



On the other hand, one finds 



|2>® 



SO in this case the expectation values agree exactly. Now we turn to 0,10,2: In the kinematical coherent 
states 



(0102)121,22) — 2:1^2 — z \z\ 
The expectation value in the physical coherent states is 



(0102) 



exp (— \zf) y/n + 1- 



\2n 



n=0 



(A.3) 



(A.4) 



The sum in this formula can not be determined in terms of elementary functions. We can, however 
study its behavior for large \z\. To this end, let us define the function 



n=l 



(A.5) 



Then we can write (|A.4]) as 



(0102) 



z\z\e-''Fs{b) 



where b = \zf. To obtain an asymptotic formula for Fa{b) for large b, we approximate the factorial 
in ([A.5|) by Stirlings formula and the discrete sum by an integral. We find 

/ ^ poo ^ 

Fa{b) ^ \ — / a;°~2 exp (— 6x(lnx — 1)) (ix 
V 2-71 Jq 
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The asymptotic behavior of this integral can be obtained by saddle point methods (see for example 
ig). We obtain 



b V 2 



(A.6) 



Therefore the expectation value ( |A.4| ) is 

{aia2)\z)^ = z\z\(^l + + 0(1^1"^)^ 



Comparing this with ( [A.3| ), we see that the expectation values in kinematical and physical coherent 
states disagree by a term of order 1. This is a small correction if 1^1 is large. 

Similar results can be obtained for more complicated functions of the Dirac observables. Consider 
for example the operator a^. It can be written as aia2A'^i~^'^^A'"2"^''^. In this case, the expectation 
value in the physical coherent states is trivial: 



Z. 



On the other hand, we find 



('^®)|21, 



Z2 



|21,22> 



e-''^e-''^F3 (61)^3(62) 



where 6,- 



Using ( [A.6| ) this can be simplified to 



Z 1 



3 1 
16 



+ 0{\z\ 



Summarizing, we find that in the simple example of two harmonic oscillators coupled by the constraint 
(|A.2|) , expectation values of Dirac observables in both, kinematical and physical coherent states, can 
be computed to any desired order of accuracy. For some observables, these expectation values agree. 
For others, there are /i-corrections. This indicates that kinematical coherent states always give the 
same answer to zeroth order as the dynamical ones and that the first corrections differ by a constant 
of proportionality of order unity so that at least qualitatively we get a good idea of which corrections 
to expect in the exact theory. 
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